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Abstract  1970 


This  two  part  paper  investigates  the  propagation  of  light  in  the  visible 
region  through  a  crystalline  medium. 

In  Pari  I,  which  is  a  slightly  shortened  version  of  the  author’s  doctoral 
thesis  of  1912,  the  objective  is  to  determine  if  the  anisotropic  arrangement 
of  ordinary  (isotropic)  dipoles  at  the  nodal  points  of  an  orthorhombic  lattice 
would  account  for  the  existence  of  double  refraction.  The  value  of  the  com¬ 
puted  birefringence  is  found  to  be  comparable  to  the  observed  value.  Addi¬ 
tionally  some  features  of  the  traditional  "theory  of  dispersion"  are  disclosed 
and  clarified.  In  the  older  theories  the  "incident  optical  wave"  plays  a  role 
but  it  is  shown  that  in  a  medium  which  extends  to  infinity  in  all  directions 
no  such  wave  should  be  assumed:  The  refractivity  is  shown  to  arise  as  an 
internal  property  of  the  medium. 

In  Part  n,  a  crystalline  medium  is  considered  as  filling  a  half-space 
and  having  a  plane  boundary  at  z  =  0.  A  plane  optical  wave  is  incident  on 
this  medium.  Because  of  the  linearity  of  the  equations  it  has  to  be  super¬ 
imposed  on  the  field  originating  in  the  crystal.  It  is  shown  that  this  inci¬ 
dent  optical  wave  is  actually  prevented  from  entering  the  crystal  because  of 
the  modification  produced  in  the  field  of  the  crystal  by  the  introduction  of  a 
boundary.  Boundary  waves  are  found  to  exist  on  both  sides  of  the  boundary. 
The  higher  their  order,  the  more  rapidly  they  attenuate  as  a  function  of  dis¬ 
tance  away  from  the  boundary  but  the  zero  order  waves  are  ordinary  un- 
dampted  plane  waves  of  vacuum  velocity  c.  The  fields  outside  and  inside 


iii 


Abstract  1970 


the  mediun-.  are  connected  by  the  Fresnel  formulae:  These  follow  from  the 
condition  that  the  optical  field  inside  the  medium  and  the  oscillations 
generated  there  form  a  self-consistent  system. 

Material  Added  1970.  The  conclusion  that  the  incident  optical  field 
cannot  penetrate  the  crystal  boundary,  together  with  a  similar  conclusion 
in  a  paper  by  Oseen,  is  the  basis  of  the  Ewald-Oseen  Extinction  Theorem. 

The  same  methods  and  field  transformations  developed  in  this  two  part 
paper  were  applied  in  two  later  papers  which  were  published  as  Parts  in 
and  IV,  These  extend  the  treatment  to  the  passage  of  X-rays  through 
crystals.  "Postscript  1970"  published  here  indicates  the  relationship  of 
these  four  parts  and  further  developments  of  the  subject  as  well  as  refer¬ 
encing  the  papers  noted  above. 


P.  P.  Ewald 


Preface 


Several  years  ago  wanting  to  knew  more  about  the  trar.'mission  of 
light  in  materials.  I  found  that  an  explanation  of  the  electromagnetic  theory 
of  transmission  of  light  in  materials  was  contained  in  the  book,  ’’Principles 
of  Optics’’  by  Bom  and  Wolf.  There  I  also  found  an  explanation  of  tt>_ 
Ewald-Oseen  Extinction  Theorem  and  reference  to  their  original  papers. 

In  order  to  have  a  better  understanding  of  the  principles,  I  attempted  to 
translate  Ewald’s  paper  which  actually  was  published  in  two  parts.  I  found 
it  to  be  not  only  basic  to  the  Ewald-Oseen  Extinction  Theorem  but  basic 
to  an  understanding  of  the  principles  of  reflection  and  refraction  of  light 
when  considered  on  an  atomic  level.  The  basic  nature  of  the  paper  is  the 
justification  for  this  translation. 

In  printing  this  translation,  the  page  format  and  the  placement  of  head¬ 
ings  has  followed  standards  established  for  AFCRL  reports.  The  numbered 
equations  of  the  original  paper  which  started  with  one  in  each  section,  have 
been  numbered  so  that  the  numbers  are  consecutive  throughout  the  transla¬ 
tion. 

The  illustrations  have  been  reproduced  photographically  from  those  in 
the  original  paper,  and  have  also  been  numbered  consecutively.  Many  of 
the  more  complicated  mathematical  equations  have  been  reproduced  photo¬ 
graphically  to  insure  against  error  in  transcription.  There  was  consider¬ 
able  use  of  "italics”  in  the  original  paper,  but  sir.ee  the  reproduction  is 
from  typed  copy  and  no  italic  characters  were  available,  italicized 
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statements  and  words  are  indicated  by  underlines.  Quantities  and  equations 
within  try.*  text  ha\e  not  been  italicized.  However,  no  confusion  should  arise 
since  the  letters  in  question  -  namely,  a,  b,  c,  i,  m,  n,  etc.,  were  not 
used  with  d  life  rent  meanings  with  normal  and  mathematical  italic  type  faces 
in  the  original  paper. 
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L.  M.  Hollingsworth 
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On  the  Foundations  of  Crystal  Optics 


introduction  to  Part  I  (Dispersion  Theory)  and 
Part  if  (Theory  of  Reflection  and  Refraction) 


The  problem  of  providing  a  foundation  of  crystal  optics  by  means  of 
the  structure  theory  of  crystals  formed  the  subject  of  my  thesis  (Ewald, 
1912^).  The  task  then  before  me  was  to  investigate  whether  dipoles, 
which,  since  the  work  of  Drude,  had  been  assumed  to  account  for  the  dis¬ 
persion  of  light,  would  also  account  for  double  refraction  when  suitably 
arranged  in  an  anisotropic  manner  so  as  to  serve  as  a  model  of  a  simple 
orthorhombic  crystal.  If  this  were  the  case,  one  would  furthermore  have 
to  determine  whether  the  magnitude  of  such  double  refraction  is  of  a  com¬ 
parable  order  to  that  observed  in  crystals  of  similar  values  of  the  axial 
ratios.  These  questions  arose,  to  be  sure,  before  the  experiments  of  Laue 
and  Ms  co-workers  on  the  interference  of  X-rays  in  crystals  had  shown 
the  correctness  of  the  structure  theory  of  crystals.  It  was,  in  fart,  hoped 
to  obtain  an  argument  in  favor  of  this  theory  from  the  theory  of  dispersion. 


^Received  for  publication  26  June  1970) 

(1)  Reference  added  1970:  Ewald  (1912)  Dispersion  und  Doppelbre c-hung 
von  Elektronengittem  (Rristallen).  (Thesis,  Univ.  of  Munich,  1912) 
G&tingen,  Dieterischscbe  Universitats  -  Buchdruckerei,  p  45. 
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The  result  of  the  thesis  is  that  without  having  to  assume  extreme  ratios 
of  the  axes,  a  lattice  arrangement  of  dipoles  can  by  itself  give  a  strong 
birefringence  because  of  the  inequality  of  the  spacing  of  the  dipoles  in  the 
different  directions.  Furthermore,  that  the  law  connecting  the  index  of 
refraction  to  the  direction  of  propagation  of  the  wave  (illustrated  by  the  sur¬ 
face  of  wave  normals)  is  fulfilled.  The  assumed  simple  model  of  the  crystal 
is  thus  correct  from  the  point  of  view  of  optics. 

The  numerically  calculated  birefringence  of  a  lattice  with  the  same 
axial  ratios  as  those  of  anhydrite  (anhydrous  calcium  sulphate),  led  to  dif¬ 
ferent  values  of  the  birefringence  than  the  observed  ones.  This  was  as¬ 
cribed  to  the  fact  that  the  crystal  consists  not  of  one  simple  crystal  lattice, 
but  of  many  simple  interpenetrating  lattices.  This  type  of  crystal  structure 
has  been  revealed  by  X-ray  investigations  and  is  now  well  established. 

It  is  not  yet  possible  today  to  decide  whether  the  birefringence  in  some 
cases  is  determined  entirely  by  the  anisotropy  of  positions,  or  whether  in 
all  cases  there  exists  also  an  anisotropy  of  the  bonding,  that  is,  a  different 
restoring  force  in  the  directions  of  each  of  the  three  principal  axes.  At 
present  we  do  not  know  of  a  birefringent  crystal  of  sufficiently  simple 
chemical  nature  for  which  the  crystal  structure  is  known  from  X-ray  dif¬ 
fraction  so  that  its  positional  birefringence  can  be  calculated  on  the  basis 
or  the  methods  used  in  my  dissertation.  Then  by  comparison  with  ex¬ 
perience,  one  could  evaluate  the  anisotropy  of  the  bonds  and  thereby  obtain 
information  about  the  internal  forces  acting  in  the  crystal.  Sulphur  might 
be  a  suitable  example  as  soon  as  its  structure  is  entirely  known.  The 
numerical  calculations  which  are  fully  given  in  the  thesis  have  been  con¬ 
siderably  shortened  in  the  present  publication. 

Our  main  interest  lay  in  the  problem  of  the  birei.\ngence:  its  magni¬ 
tude  and  the  way  it  comes  about.  However  it  turned  out  that  some  novel 
facets  appeared  from  the  point  of  view  of  the  general  theory  of  dispersion. 
The  precisely  given  positions  of  the  dipoles  made  a  more  rigorous  calcula¬ 
tion  possible  than  in  past  treatments  of  the  dispersion  of  isotropic  bodies 
which  required  averaging  procedures.  Small  conceptual  differences  here 
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forced  attention  upon  ti.  .mselves.  The  first  questions  is:  How,  exactly, 
are  the  outgoing  fields  of  the  individual  dipoles  to  be  summed  in  the  crystal 
which  we  think  of  as  being  of  infinite  size  ?  This  is  a  question  of  conver¬ 
gence  which  however  remains  meaningful  also  for  the  finite  crystal. 

A  second  new  aspect,  at  least  with  respect  to  the  Planck  and  Lorentz 
theory,  was  the  omission  of  the  "incident  wave".  Planck  and  Lorentz  let 
the  single  dipole  be  excited  by  an  electrical  force,  which  consists  of  two 
parts,  the  'External  one"  (G0)  and  a  second  contribution  which  comes  from 
the  other  dipoles  of  the  body.  The  first  part,  ©j,  ,  found  no  place  in  the 
theory  of  the  infinite  crystal.  The  mutual  interactions  alone  maintained 
the  dipole  vibrations,  so  that  the  complete  state  of  oscillation  of  the  crystal 
had  to  be  looked  upon  as  3ome  type  of  free-vibration,  not  requiring  an 
external  forcing  field.  The  idealization  of  the  problem  by  considering  a 
crystal  with  no  boundary  was  found  to  contain  justification  for  the  omission 
of  the  exterior  exciting  field  ©0  .  In  a  finite  crystal,  its  boundary  should 
play  the  role  of  deleting  the  incident  wave  ©q  for  the  entire  interior  of  the 
crystal.  Only  in  this  way  was  it  possible  to  understand  that  in  an  unbounded 
crystal  no  exterior  excitation  can  be  noticed. 

In  this  assertion  lies  the  admission  that  for  a  complete  understanding 
of  dispersion  theory  the  treatment  of  reflection  and  refraction  is  absolutely 
essential.  This  supplementary  theory  of  the  events  in  the  boundary  of  th*- 
body  will  be  given  in  the  second  part  of  this  paper,  and  that  part  is  new 
and  original.  '  It  treats  the  following  problem: 

Given  a  bounded  crystal  -  which  fills  the  entire  upper  half  space  above 
the  boundary  surface  z  =  0;  from  below  a  plane  wave  ©0  is  incident  on  the 
boundary  at  an  arbitrary  angle.  Find  the  mode  of  vibration  of  the  dipoles, 
and  whether  this  leads  to  the  formation  of  refracted  and  reflected  waves  - 
and  implicitly  to  the  Fresnel  formulas.  Finally,  can  the  anticipated  in¬ 
fluence  of  the  boundary  be  shown? 

(I)  Note  added  on  galley  proof.  Mr.  C.W.  Oseen  recently,  in  a  dif¬ 
ferent  context  has  also  called  attention  to  the  necessity  to  treat  dispersion 
theory  in  the  sense  discussed  above.  Oseen,  C.W.  (1915)  Physikal. 
Zeitschr.  16:404. 
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This  problem  of  the  bounded  crystal  appears  at  first  sight  to  be  con¬ 
siderably  more  difficult  than  that  of  the  unbounded  crystal.  In  actual  fact, 
however,  the  vibrational  mode  as  found  in  the  thesis  which  resembles  the 
free  vibrations  of  a  mechanical  system  is  a  stepping  stone  for  dealing  with 
the  bounded  crystal.  Just  because  the  surface  layer  screens  the  inside  of 
the  crystal  from  the  invasion  of  the  external  wave,  a  mode  of  oscillation 
inside  the  crystal  becomes  possible  as  if  there  were  neither  a  boundary  nor 
an  incident  wave.  Since  this  mode  of  vibration  is  well  known  to  us  for  the 
unbounded  crystal,  it  requires  little  effort  to  reduce  the  case  of  the 
bounded  crystal  to  the  case  of  the  unbounded  crystal. 

The  second  part  of  this  paper  may  be  considered  an  extension  or  con¬ 
tinuation  of  the  problem  treated  in  the  thesis  and  in  particular  as  a  proof 
of  a  statement,  made  there.  Furthermore,  a  particular  point  that  was  not 
treated  before  will  be  discussed.  This  is  the  transition  from  the  precise 
microscopic  field  that  exists  between  the  dipoles,  to  the  macroscopic 
average  field  that  is  used  in  the  phenomenological  theory  for  the  treatment 
of  reflection  and  refraction  and  whose  amplitudes  are  determined  by  the 
Fresnel  formulas.  These  are  concerned  with  the  comparison  of  the  field 
strengths  outside  and  inside  of  the  crystal.  One  has  therefore  to  watch 
closely  the  relations  between  the  calculated  microscopic  field  in  free  space 
between  the  oipoies  of  the  body  and  the  Maxwellian  field  quantities  £  and 
X  inside  the  body.  For  the  dispersion  problem  proper,  we  are  only  in¬ 
terested  in  establishing  the  phase  velocities  with  which  the  quantities  £ , 

%  and  propagate.  For  this  purpose  it  is  not  necessary  to  distinguish 
between  6  and  3) ,  and  according  to  the  thesis,  the  results  can  be  obtained 
solely  from  the  vector  potential  without  even  specifying  the  field  strengths. 
This  discussion,  however  becomes  necessary  in  connection  with  the  Fresnel 
formulas.  Section  5,  Part  n,  really  belongs  just  as  much  to  Part  1  (dis¬ 
persion)  as  to  the  theory  of  refraction,  but  was  omitted  from  Part  I  in 
order  not  to  interrupt  the  train  of  thought  for  the  treatment  of  dispersion 
theory. 


o 


Finally,  let  it  be  mentioned  that  the  behavior  of  crystals  for  very  short 
waves  (x-rays) ,  is  contained  in  the  rigorous  calculations.  See  Section  7. 

It  is  of  particular  interest  to  follow  the  transition  from  the  case  of  x-rays 
with  a  multiplicity  of  interference  rays  to  the  case  of  optical  rays  with 
only  three  waves  (the  incident,  the  refracted,  and  the  reflected  waves). 
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Part  I.  Dispersion  Theory* 


1.  STATEMENT  OF  PROBLEM 

(1)  Given  a  simple  orthorhombic  If  i.Ice  with  lattice  points  at  cartesian 
coordinates 

X  —  2al ;  Y  —  2bm  ;  Z  =  2cn 
(l,  m,  n  integers  from  ~  CD  to  +  oc) 

Electrons  are  bound  to  these  lattice  points  with  quasi-ela  tic  and  isotropic 
forces.  The  last  statement  means  that  the  restoring  force  between  the 
electron  and  the  lattice  point  depends  only  on  the  magnitude  but  not  on  the 
direction  of  the  displacement.  With  what  velocity  and  what  attenuation  can 
one  have  a  plane  wave  of  a  given  frequency  propagate  in  this  medium  ? 

(2)  We  begin  with  the  reminder  that  for  the  purpose  of  calculating 
optical  fields  the  oscillating  electrons  of  charge  -e  can  be  replaced  by 
oscillating  dipoles.  We  imagine  fixed  charges  +e  and  -e  to  be  placed  at 
the  lattice  points.  Then  +e  together  with  the  displaced  electron  form  the 


*  Shortened  thesis  Munich  1912.  I  wish  to  not  fail  to  express  rry 
sincere  thanks  in  this  place  to  Profes  sor  Sommerfeld  for  his  suggestions 
and  support  in  this  work. 
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dipole,  while  the  remaining  charges  -e  at  the  lattice  points  produce  an 
electrostatic  field  in  which  we  are  not  interested. 

In  the  equation  of  a  plane  wave: 

6  «  WeiEe-"'*  •'**), 

E  is  the  amplitude,  n  the  frequency,  and  s  is  the  coordinate  measured  in 
the  direction  of  the  wave  noimal.  The  sought  for  values  of  phase  velocity 
and  attenuation  are  contained  in  the  constant  k  whose  real  part  gives  the 
phase  velocity  q  =  n/k^,  while  the  imaginary  part  kg  produces  the  attenua¬ 
tion. 

One  now  makes  the  assumption  that  the  planes  perpendicular  to  the 
planes  s  =  constant  are  planes  of  equal  dipole  phase.  The  phase  velocity  q 
with  which  the  phase  propagates  in  the  lattice  is  at  first  not  fixed.  One  now 
shows  that  the  superposition  of  the  outgoing  individual  fields  of  the  dipoles 
pro'iices  a  (continuous)  field  strength  which  has  likewise  the  form  of  a 
plane  wave;  this  being  a  direct  consequence  of  the  assumption  made  for  the 
(discrete)  dipole  oscillations:  furthermore,  that  the  frequency  and  attenua¬ 
tion  of  both  waves  -  of  the  field  and  of  the  dipole  excitation  -  are  identical. 
In  order  to  determine  these  quantities  which  are  both  contained  in  the  com¬ 
plex  constant  k,  we  have  to  insure  the  dynamical  possibility  or  consistency 
of  the  original  assumption.  This  requires  that  each  dipole  is  excited  into 
oscillation  by  the  field  impinging  on  it.  This  condition  can  be  fuililleu  for 
any  frequency  by  proper  choice  of  the  arbitrary  constant  k  and  this  yields 
the  relationship  between  frequency  and  phase  velocity,  namely  the  disper¬ 
sion  formula. 

2.  METHOD  OF  CALCULATION 

Let  XYZ  be  the  coordinates  of  a  dipole  and  xyz  the  coordinates  of  a 
field  point;  and  further  let 


S  =  X  cos  (*,  x)  +  V  cos  (s,  y)  +Z  cos  ($,  z) 
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be  the  distance  of  a  dipole  from  the  plane  of  equal  dipole  phase  which 
contains  the  origin:  let  q  La  the  (undetermined)  dipole  phase  velocity,  n 
the  frequency  in  the  time  factor  c~  1  so  that 

e  ~  *'*'  «• 5 


is  the  phase  factor  propagating  in  tie  direction  -S  of  the  dipole  excitation. 

The  individual  field  emitted  by  a  dipole  with  the  charge  =e  is  then 
characterized  by  the  Hertz  vector 


ea’i 


(i) 


in  which  R  =  ((X-x)2  +  (Y-y)2  +  (Z -z)2]  1^2,  signifies  the  distance  from 
dipole  to  field  point.  The  term  R/c  in  the  exponent  indicates  that  the  field 

itself  propagates  toward  the  field  point  with  the  same  velocity  c  as  in  free 

* 

space. 

Since  all  dipoles  have  the  factor  e  1  in  common,  we  omit  it  and  in¬ 
troduce 


-  -  -  -  i 


(2) 


and  then  write  for  the  significant  part  of  the  Hertz  vector 


»  _  1  iKS-ikS 

P  *»  ea  •  £  •  e  , 


(3) 


and  recall  that  in  a  material  with  dielectric  constant  1,  the  electric  field 
strength  6  and  the  magnetic  field  strength  £  can  be  obtained  through  the 
operations 


♦The  original  reads  ’)pure  ether"  but  free  space  has  been  substituted 
as  the  modern  equivalent  expression. 
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coming  from  the  test  dipole.  The  equation  for  the  oscillation  then  is  given 
[from  Eq.  (4)j 

nr  e*curl  curl  (a  Tl'  *  (S) 

Here 

p  =  a  • 

is  the  elongation,  m  is  the  mass,  e  is  the  charge,  f  is  the  quasi-elastic 
force,  g  the  coefficient  of  friction.  By  introducing  the  proper  frequency 
of  the  dipole 


the  equation  becomes 

a  n1  4-  *0*  —  in  *  j  =  ^curl  curio  FT. 
To  shorten  this  we  set 


and  replace  curl  curl  with  its  equivalent:  (grad  div  -A).  Thus  we  have 
to  fulfill 

a  (Q  ~  All')  =  grad  div  q  77'  (10) 

or  in  components  (the  indices  cm  the  potential  U*  indicate  differentiation 
while  those  on  the  amplitude  a  signify  component  formation): 
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r,?2  +  j  n  -  n\ 3  -  Qr  -  o,  nst  =  i» 

-  o,  /?•,,  +  Q,{--  +  J  77'  -  Z7',,)  -  q,  /7  -  0 

-  a, /7xl  -  ar  Tl'it  +a,(2+ J/7'-  77'. J  r,  u .  (10a) 

The  condition  for  simultaneous  fulfillment  of  these  three  homogeneous  equa¬ 
tions  is: 


\Q  +  AW  -  n\ .  -  n\ 


-  n\ 


i  —  / 1 

-  n\ 


ii  4.  J  /7'  -  ntj  -  I2'ft  i  =  U. 

-  77*  a  +  Air-  //'„! 


The  13'. . .  depend  on  the  constant  k  which  determines  the  dipole  phase 
and  this  is  the  same  as  that  which  determines  the  propagation  of  the  optics! 
field.  Thus  the  determinant  Eq.  (11)  expresses  the  relationship  between 
the  index  of  refraction 

r  t  _ 

'"7 -X  m 

and  the  quantity  Q,  which  contains  the  electron  constants  and  the  frequency 
of  the  optical  wave.  *Law  of  dispersion.) 

We  obtain  the  principal  indices  of  refraction  of  the  orthorhombic  lat¬ 
tice  by  assinning  the  direction  of  the  dipole  oscillation  to  be  along  one  of 
the  three  axes;  for  each  of  these  axes  this  can  be  done  with  two  different 
directions  of  propagation  of  the  optical  wave.  We  thus  obtain  in  the  first 
Instance  six  principal  Indices  of  refraction,  which  depend  on  the  condition 
of  polarization  and  on  the  direction  s  of  the  wave  normal: 

Oscillation  a=a„  Propagation  along  z  :  k„,  rm 
•  -  •  b':  km 


Q  —  Q, 


Q  —  0, 


X . 

*  • 

^ 

y-h*  v 


Herein  lies  a  departure  of  the  lattice  from  the  actual  crystal  in  which  oniv 
three  principal  indices  of  refraction  are  observed,  and  we  will  now  con¬ 
sider  the  origin  of  the  six  indices  in  a  more  direct  way. 

Let  the  dipole  amplirude  be  in  the  direction  of  the  x-axis.  This  leaves 
us  two  different  ways  to  dispose  of  the  plane  of  constant  dipole  phase  (see 
Figure  1): 

(1)  The  phase  plane  is  the  xy-plane  (propagation  along  z). 

(2)  The  phase  plane  is  the  xz-plane  (propagation  along  y).  Or:  account 
of  the  non-equivalence  of  the  y-  and  z-  directions  in  an  orthorhombic  lattice, 
these  two  cases  lead  to  two  difftrent  exciting  fields  and  there  are  in  general 
actually  six  prip.cipsl  ipdicss  of  rsfrsction  to  distinguish* 

It  is  clear,  however,  without  further  calculation  that  the  difference 
between  the  two  indices  of  refraction  v  and  »-  is  very  slight.  For  it 
arises  from  the  fact  that  a  phase  difference  exists  between  the  test  dipole 
and  its  neighbors  -  in  one  case  the  neighbor  in  the  z-direction,  in  the  other 
case  those  along  y.  The  phase  difference  between  adjacent  dipoles  being 
very  small  anyway7,  because  of  the  large  wavelength  of  the  optical  field,  any 
effect  produced  by  the  inequality7  of  the  two  phase  differences  will  be  ex¬ 
tremely  small. 

ft  will  therefore  be  sufficient  henceforth  to  distinguish  only  the  three 
principal  refractive  indices  t*v,  v  which  refer  to  the  directions  of 
oscillation  a,,  a,,  a,  respectively.  Their  differences  determine  the  bire¬ 
fringence.  The  differences  between  them  are  considerably  greater  than 
those  caused  by  the  directions  of  wave  p*-opagation.  because  thev  are  nearly 


Figure  1 
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independent  of  the  wavelength.  Figure  la  shows  the  dipoles  with  their  os¬ 
cillations  in  the  x-direction  and  Figure  lb  the  same  lattice  with  the  oscil¬ 
lations  in  the  y-direction.  For  an  understanding  of  the  physical  significance 

of  the  difference  between  v  and  v  .  (optical  rotation  power) ,  see  the  con- 

xy  xz 

eluding  portion  of  Section  4. 


When  computing  the  principal  indices  of  refraction  (propagation  along 

one  axis  and  oscillation  along  another  axis),  the  three  coordinate  planes 

of  the  lattice  remain  planes  of  symmetry  and  for  this  reason  at  the  lattice 

points  the  differential  quotients  nf  =  II'  =11*  -  0;  and  thus  the  deter- 

xy  xx  yz 

minant  reduces  to  the  diagonal  terms  and  one  obtains  from  Eq.  (10a)  the 
equation  for  the  determination  of  v 

XZ 


t h  4“  {All'  —  /7'j,),  =  0. 


(13) 


The  index  z  is  to  remind  one  that  the  specialized  formulas  for  propagation 
along  z  must  be  used  in  order  to  find  v  .  .  In  this  equation  the  second 

XZ 

derivatives  of  the  potential  -  as  the  approximations  show  -  contain  only  the 
ratio  k/kQ  =  i\  but  not  the  frequency  n  which  appears  only  in  the  electron 
constant  fl. 

If  one  eliminates  the  electron  constant  O,  from  Eq.  (13)  and  the  similar 
equation  for  v 

Q  +  (AW  -  W„)t  =  0  (13a) 


one  obtains  an  implicit  relationship  between  the  two  principal  indices  of 

rc,>  ction  v _ and  v  which  does  not  depend  on  the  frequency  n;  nor  does 

- xz  yz - 

it  depend  on  the  number  and  characteristic  constants  of  the  dipoles;  this 
relationship 


(AW  =  (AW  -  W„\. 


(14) 
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depends  solely  on  the  structure.  In  the  case  of  an  orthorhombic  lattice 
this  equation  leads  to  a  relationship  previously  found  by  T.  H.  Havelock 
(see  Section  4) : 


Const. 


whose  non-dependence  on  the  wavelength  Havelock  found  to  be  confirmed 

in  some  crystals.  The  calculation  of  11*  and  II'  enables  us  to  determine 
-  xx  yy -  — 

the  magnitude  of  the  measure  D  of  the  structural  birefringence  for  a 
given  axial  ratio. 


3.  SOME  GENERAL  REMARKS  ON  PROCEDURE 

(1)  The  summation  of  the  individual  potentials 


(15) 


which  is  required  in  the  previous  section  for  the  calculation  of  the  optical 
field  as  well  as  of  the  "exciting  field, "  has  no  definite  meaning  for  real 
values  kg  =  n/c  and  k  =  n/q.  Because  of  the  oscillating  factors  in  the 
numerator,  the  decrease  of  the  partial  terms  as  l/R  is  just  about  sufficient 
to  make  possible  a  finite  value  of  the  sum;  however,  this  value  is  indefinite 
and  dependent  upon  the  manner  of  the  execution  of  the  summation  because 
the  convergence  of  the  series  is  "conditional".  To  arrive  at  an  unequivocal 
result,  the  observations  of  the  previous  sections  have  to  be  supplemented 
by  summation  instructions  which  have  to  be  derived  from  the  nature  of  the 
physical  problems. 

To  start  with,  one  is  relieved  of  all  difficulties  of  convergence  by  con¬ 
sidering  a  finite  piece  of  crystal  which  consists  of  a  large  number  of  di¬ 
poles  N  .  When  a  light  wave  falls  onto  this  crystal,  the  dipoles  begin  to 
oscillate  and  create  together  with  the  incoming  wave  (£&  ,  an  electric  field 
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in  the  body.  Tins,  in  turn,  determines  the  manner  of  the  oscillations,  the 
propagation  velocity,  etc.  .  as  shown  in  the  preceding  sections. 

In  the  finite  body,  this  field  depends  essentially  on  the  boundaries.  If 
the  surface  on  which  the  incidence  occurs  is  plane,  the  wave  is  refracted. 

If  it  is  undulating,  diffraction  spectra  of  different  orders  develop  and  all  of 
the  diffraction  phenomena  at  the  boundaries  of  the  body  are  within  the  solu¬ 
tion  of  this  problem.  The  dependence  on  the  surface  extends  throughout  the 
entire  inside  of  the  body. 

Obviously,  the  meaning  of  the  prrblem  of  dispersion  is  more  restricted 
than  the  proposed  task:  it  is  enough  to  determine  the  propagation  velocity  of 
a  wave,  that  is,  a  material  constant  which  is  independent  of  outside  condi¬ 
tions  (Doundaries) .  The  very  existence  of  this  problem  shows  that  the 
essential  process  of  wave  propagation  at  a  point  cf  tie  body  is  caused  by 
the  environment  of  this  point.  Neither  the  boundaries  nor  the  outer  wave 
can  influence  this  process  -  or  else  it  would  have  to  vary  in  the  body 
from  place  to  place.  The  progress  of  a  wave  Inside  of  a  body  is  a  dynami¬ 
cally  closed  process  which  inherently  has  the  capacity  that  when  excited  in 
one  part  of  the  body,  it  can  move  on  to  other  daces  without  the  help  of  ex¬ 
ternal  forces.  The  variety  of  the  diffraction  phenomena  of  a  bounded  body 
results  because  the  boundary  initiates  many  types  of  oscillations  in  the 
body.  These  then  continue  to  develop  in  a  ’’dynamically  closed"  manner. 

If  our  aim  is,  in  the  pure  theory  of  dispersion,  to  find  the  'dynamically 
closed"  state  of  oscillation,  then  we  have  to  eliminate  the  influence  of  the 
boundary.  We  must,  because  a  material  constant  is  involved,  consider 
only  the  influence  of  the  environment  at  the  place  of  the  propagation.  That 
the  environment  should  extend  everywhere  over  many  wavelengths  does  not 
enter  into  che  question  of  convergence.  The  simplest  way  to  proceed  is  to 
return  to  the  infinite  crystal,  but  to  introduce  instead  of  the  everywhere 
equally  strong  dipole  amplitude  a  a  gradually  diminishing  amplitude  at 
distances  fer  from  the  test  point 


aa~ 
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If  x^'l  /  A.  this  produces  no  change  within  many  wavelengths  from  the 
test  point.  However,  it  assures  the  determinancy  of  the  problem  and  the 
unconditional  convergence  of  the  sum  of  the  individual  potentials.  This 
sum  may  now  be  written 

-xti  ik^lt-  ikS  _  ik/it-iks 

2“  i  -2ai 

and  is  formally  identical  with  the  old  sum  except  for  the  fact  that 
hg  —  Atq  H"  ix 

has  the  imaginary  part  *  . 

In  the  case  of  an  absorbing  crystal,  for  which  one  would  have  to  as¬ 
sume  an  amplitude  growing  exponentially  in  the  direction  opposite  to  that 
of  propagat^n  (k  =  +  ikg)  definite  values  of  the  sums  can  evidently  only 

be  obtained  by  assuming  a  value  of  sr  larger  than  For  x  -  k2  would 
Just  compensate  for  the  increase  of  dipole  amplitude  and  reduce  the  sums 
to  the  same  indefiniteness  as  in  the  non-absorbing  crystal  and  x  -  0  . 

In  the  following  we  re-introduce  m  place  of  kg  explicitly 

&io  ==  -f-  * 

whereby  k^9{=x)  >  .  The  result  of  the  summation  with  this  complex  kg, 

is  one  of  the  possible  sum  values  with  real  kg.  By  executing  the  summation 
with  a  complex  value  and  reducing  its  imaginary  part  to  zerc  in  the  sum,  a 
value  of  the  sum  is  obtained  which  is  free  of  any  influence  of  the  boundary. 

(2)  From  the  last  statements,  it  sliould  not  be  surprising  that  in  the 
foregoing  paragraphs  no  mention  was  mace  of  any  wavs  (50  incident  from 
the  outside.  The  force  which  stands  on  the  righthand  side  of  the  equation 
for  the  oscillations  contains  only  fields  radiated  by  the  other  dipoles  of  the 
body,  not  as  with  Planck  and  H.A.  Lorentz,  a  term  cov.-esponding  to  an 
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incident  wave.  Such  a  term  would  contradict  the  interpretation  of  the  wave 
propagation  as  a  "dynamically  self-contained"  phenomenon.  The  impos¬ 
sibility  of  the  existence  of  the  primary  wave  inside  the  body  is  evident 
from  the  result  of  the  summation  of  the  individual  fields:  the  assumed 
dipole  oscillations  produce  a  field,  which  by  itself  has  the  form  of  a  plane 
wave  travelling  with  velocity  q.  If  this  field  were  to  be  superimposed  on 
the  primary  field  %  which  has  another  direction  and  travels  with  another 
velocity  c,  these  fields  would  never  combine  within  the  body  to  form  a 
single  plane  wave. 

According  to  the  point  of  view  in  subparagraph  (1),  it  is  dear  that  the 
absence  of  Go  must  be  due  to  an  effect  of  the  boundary.  The  plane  boundary 
of  a  body  fulfills  the  functions  of 

a.  compensating  the  incident  wave  everywhere  in  the  interior  of  the 
body; 

b.  initiating  the  "dynamically  self-contained"  mode  propagation  of 
the  refracted  wave; 

c.  producing  the  reflected  wave  in  outer  space. 

A  forceful  proof  for  the  required  effect  of  the  boundary  may  be  obtained 
from  considering  a  ray  of  limited  cross  section  falling  under  an  angle  on  a 
plane  boundary  (Figure  2).  But  for  the  fact  that  the  exposed  part  of  the 
boundary  produces  a  plane  wave  of  equal  direction  and  phase  velocity,  but 
opposite  phase,  as  the  incident  wave,  there  would  have  to  exist  an  optical 
field  also  in  the  region  between  the  dashed  lines  of  the  figure. 


Figure  2 
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4.  RELATIONSHIP  TO  OTHER  THEOR'ES 

This  work  comes  closest  to  the  dispersion  theory  of  Lorentz  and 
Planck.  The  main  difference  is  that  in  their  treatment  of  the  isotropic 
body,  average  values  are  used  whereas  in  this  paper  strict  consideration 
of  the  geometric  location  of  the  dipoles  is  employed. 

One  further  difference  of  the  basic  assumptions  has  alreaay  been  men¬ 
tioned  above;  namely,  that  we  omit  tlie  incident  wave  .  It  remains  to 
say,  however,  that  the  results  of  the  above  mentioned  theories  are  con¬ 
firmed  in  this  paper  when  the  results  are  specialized  for  a  cubic  crystal. 

In  computing  the  birefringence,  this  work  exceeds  the  aims  c  f  the  above 
mentioned  theories  and  is  closely  connected  with  a  group  of  papers  by 
Lord  Rayleigh^,  T.  Havelock^,  and  N.  Kasterin^.  These  papers  start 
out  from  Rayleigh’s  investigation  of  the  properties  of  a  medium  in  which 
obstacles  in  the  form  of  small  spheres  of  different  dielectric  constant  are 
imbedded  in  a  lattice  arrangements.  Only  the  static  problem  of  finding  the 
potential  distribution  is  treated.  On  this  basis  Havelock  uses  the  dielectric 
constant  for  the  composite  medium  and  assumes  that  it  varies  with  wave¬ 
length  in  accordance  with  Cauchy’s  formula.  This  gives  him  a  dispersion 
formula. 

Kasterin,  however,  has  succeeded  with  a  method  similar  to  those 
used  by  Rayleigh  to  find  a  solution  for  the  wave  equation  Au  4-  k*u  =  0 
which  fulfills  the  boundary  condition  on  the  surfaces  of  the  spheres  and 
which  represents  plane  waves.  It  is  a  pity  that  this  excellent  piece  of  work 
has  only  been  applied  to  acoustics  problems. 

The  connection  between  Kasterin’s  work  and  the  present  study  is  based 
(4) 

upon  the  remark'  ,  that  the  way  a  plane  wave  is  affected  by  a  small 

(1)  Lord  Rayleigh  (1892)  Phil.  Mag.  34. 

(2)  Havelock,  T.  H.  (1906)  Proc.  Roy.  Soc.  London  77:170. 

(3)  N.  Kasterin,  Lorentz-Jubelband,  or  Versl.  Acad.  Wetensch. 

Amsterdam  460:1897-1898. 

— . . — 

(4)  See  P.  Debye,  Der  Lichtdruck  auf  Kugeln,  diss.  p.  50. 


spherical  obstacle  can  be  replaced  in  the  first  approximation  by  super¬ 
position  of  a  dipole  spherical  wave  issuing  from  the  obstacle. 

in  a  second  paper,  Havelock^  has  approached  birefringence  with  an 
idea  occurring  in  Planck’s  theory,  namely  that  in  isotropic  bodies  there 
exist  3  around  each  test  dipole  a  spherical  cavity  such  that  the  molecules 
lying  within  it  do  not  on  the  average  excite  the  test  dipole  in  the  middle  of 
the  cavity.  According  to  Havelock,  there  exists  in  a  crystal  an  ellipsoid 
of  similar  properties.  The  birefringence  can  then  be  expressed  by  its 
eccentricity,  but  it  is  not  connected  to  measurable  quantities.  Havelock 
finds  in  his  work,  the  relationship  D„  (Section  2),  between  the  principal 

refractive  indices  which  should  be  independent  of  the  wavelength. 

(2) 

Langevin  comes  to  the  same  result  by  assuming  that  each  molecule 
is  anisotropic.  Thus,  if  it  is  through!  of  as  being  isolated  in  space,  it  is 
assumed  to  have  three  different  P’gen-frequencies.  We  shall  show  in 
Section  9  that  the  qualitative  result  of  the  assumption  of  an  anisotropic 
molecule  must  be  the  same  as  that  resulting  from  the  assumption  of  a  lat¬ 
tice  structure. 

NOTE:  Since  the  first  publication  of  this  paper  (thesis),  the  same  problem 
has  also  been  treated  by  Max  Bern  in  his  "Dynamik  der  Kristallgitter" 
(Teubner,  1915).  Bom’s  mathematical  method  follows  closely  the  presen¬ 
tation  which  Hilbert  gave  of  the  same  subject  in  his  lectures  on  optics  in 
1913.  In  these  a  consistent  mathematical  development  of  potentials  ac¬ 
cording  to  powers  of  the  small  quantities  1,  a,  p.y  etc. ,  is  carried  out. 
Bora  has  treated  the  problem  in  a  very  general  form,  namely  for  non- 
orthogonal  lattices  containing  different  sorts  of  electrons  or  ions.  He  also 
assumes  that  the  charges  are  not  bound  elastically  to  fixed  lattice  points, 
but  only  to  one  another,  and  that  therefore  mechanical  oscillations  occur 
simultaneously  with  electrical  ones.  The  case  treated  in  my  thesis  can  be 
derived  from  Bora's  equations  by  assuming  an  orthorhombic  lattice  in 
which  all  simple  lattices  from  which  it  is  built  have  infinite  mass  except 
one.  The  one  remaining  type  of  charge  oscillates  then  as  though  it  were 
bound  to  fixed  lattice  points. 

(1)  Havelock,  T.  (1907)  Proc.  Roy.  Soc.  London  80:28. 

(2)  Langevin,  P.  (1910)  Sur  1  is  birefringences  electrique  et 
magnetique,  Le  Radium  7, 
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converges  wben  £  is  so  behaved  that  for  its  small  values  <0,  and 

\*/ 

for  large  values  Sc  (a*)  >  0  .  This  is,  for  example,  the  case  when  €  runs 
out  from  the  origin  within  the  sector  -  r/4  +  9  and  -  */4  and  continues  in  a 
straight  line  to  infinity. 

The  value  of  the  integral  may  easily  be  determined  (see  Riernann- 
Weber,  d.  part.  Dgl.  a.  math.  Pfaysik  L,  5th  Edition,  p.  23) : 


-»*  +  p  r/^ 

*de  =  ¥£e*2il 


(17) 


By  putting  e  -  R-  n,  the  path  of  integration  remains  essentially  unchanged 
and  we  have 


r  -*•*»+ 

J‘ 


*  l 


t  s 


Finally  by  writing  (l/2ik  R  instead  of  k,  we  obtain  the  single  integral 


He, 
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whose  path  can  be  deformed  sc  as  to  coincide  with  the  positive  real  axis, 
except  near  the  origin,  where  it  has  to  lie  within  the  angle  6.  ( See  Figure  3'. 


Figure  3 


5.2  The  Triple  integral 


In  order  to  obtain  a  representation  in  which  the  exponent  of  e  is  a 
linear  function  of  x,  y,  z,  the  simplest  wav  is  to  transform  the  single 
integral  representation  Eq.  (19)  as  follows: 

We  have 


*r" 

t  dl-t-*'’ J  , 


2ee-^‘  V'* 


We  substitute  in  Eq.  (IS)  in  place  of  e  ~  r*’  the  integral  over  ^  standing 
on  the  left  hand  side  of  this  equation  and  similarly  for  c~r'r'  and  «--■*** 
analogous  integrals  over  the  variables  p  and  v.  Thus  we  obtain  the 
quadruple  integral: 


>{1*4-  +  »*>  +  — 


(20) 


The  order  of  integrations  may  be  changed  and  the  integration  over  r 
performed  first: 

ft  «-»<*« 

m 

This  leaves  the  triple  integral  standing 


(22) 


which  integral  has  the  desired  very  suitable  form  for  summation. 
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II  is  seen  that  the  substitution  -  \  . .  in  place  of  +  \  . .  changes  nothing 
but  the  sign  of  the  exponent.  Furthermore,  if  R  is  the  distance  of  the  point 
xyz  from  the  point  XYZ,  rather  than  from  the  origin,  the  generalized  in¬ 
tegral  is  valid 


di  du  dr 


(23) 


Finally  it  should  be  noted  that  it  is  possible  to  shift  the  integration  path 
from  the  real  axes  provided  that  tie  condition 

k*  -  «  -  p*  -  r*  =0 

is  avoided.  In  particular,  if  the  positive  imaginary  part  of  k  is  nonzero, 
identical  integration  paths  a  snort  distance  above  the  real  axis  may  be  em¬ 
ployed  for  the  three  variables.  (That  is,  between  the  axis  and  the  point  k). 


6.  THE  TOTAL  POTENTIAL  AND  THE  OPTICAL  FIELD 

In  order  to  obtain  the  total  potential  of  the  ortho  bom  bie  point  lattice, 
we  first  make  use  of  the  triple  integral  [Section  5,  Eq.  (23)]  and  sum 
the  individual  potentials 


, — -i*5 


_«»  r  r  (•  ±:i(Z-4)±lr(T-y',±i*iZ-z) 


(24) 


The  coordinates  of  the  lattice  points  are: 


X=  '2*1  Y  —  2  bitu,  Z  =  '2cn, 


where  I,  m,  n  take  on  all  integral  positive  and  negative  values,  (a,  b,  c) 
are  the  half  distances  between  neighboring  points  on  the  crystal  axes  and. 
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at  the  same  time,  the  spacings  of  adjacent  planes  of  symmetry  in  the 
crystal,  is  the  propagation  constant  of  the  individual  spherical  wave, 
k  that  of  the  dipole  phase.  We  write 


kS  =  k(X  cos(S,  x )  +  lr  cos(S,  y)  +  Z  cos (S,  z)) 
—  aX  +  P  Y  4-  yZ ; 


(25) 


whereby 

a*  -f  3X  -*■  y*  =  kz.  (25a) 

The  total  potential  then  is 


— «o . .  +«o 


V  -  ~  ~ 


(26) 


We  carry  out  the  summation  under  the  integral.  The  triple  sum  is  then  the 
product  of  three  linear  geometric  series  of  winch  the  first  is 


Convergence  of  this  sum  for  positive  and  negative  values  of  £  can  be 
achieved  thanks  to  the  uncertainly  of  the  sign  of  \  provided  that  the  im¬ 
aginary  part  kgg  is  greater  than  kg  (see  Section  3).  For  if  this  is  so,  we 
can  let  X  rsm  cm  a  line  parallel  to  the  real  axis  and  above  <?,  similarly  for 
n  and  v  with  imaginary  parts  greater  than  those  of  8  and  y*  respectively. 
Thus  +X  -  a  has  a  positive  imaginary  part  which  insures  that  the  sum  con¬ 
verges  for  positive  values  of  i.  The  sum  over  negative  l  converges  when 
we  reverse  the  sign  of  X . 


*  One  finds  the  proof  in  the  thesis  note  1. 
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The  summation  over  t  therefore  yielJs 


T  3 


—  oo 
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1  +  a*1**  —  S  cos  2  aa  a3,al 
#  —  U(*-rl«)  _  ,»Aa  +  2laa  _  4U(r  +  2a) 

2 co«2 a i  -  2eo*  2  aa  " 

_  ■l‘n^t  +  go)  -  •2,mm  cn  l*  ' 

~  ec^2ai-oo*2aa 


(27) 


Appropriately  performing  the  summations  over  m  anu  n,  the  Total  Potential 
thus  is 


i  rrr»inl(z+*o)-s3im*ilnlx  »ittfi(y-f>fr)-*3i<>*aa»»y 
©o*  2«i  —  eo*  *•«  ’  CO*  Tbfi  —  00* 

■jayOH-l*,  —  *3*  ***!■»«  .  i *  dft  4v  f 

catiye  —  ooalej  kj  - 1* -  p*  -  ** 


where  the  integration  paths  for  (A  flu)  are  lines  parallel  to  the  real  axes, 
and  lying  above  {apy),  respectively.  (Path  1,  Figure  4). 
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Figure  4 


Let  us  remove  a  part  of  the  integrand.  We  write  in  the  A-integral 


riaKs-f  2o)-«2<**»inXa? 
co* 2 a i  -  cmS ay 

*inl:rco*?al  +  co*  l  x  »in  8  a  l  -  *au,*lals 

m  co*  2a  i  -  co*  2a  a 


—  sin  A*  *f 


co»  l  x  »in  2a  i  —  t  »in  X x  ito  3  a  a 
co*  2a  1  -  co*  2a a 


Here  the  first  term,  sin  Ax,  has  lost  the  denominator  and  with  it  a  pole  of 

the  integrand.  For  this  part  of  the  integral  we  can  then  shift  the  path  of  A 

2  2  2  2 

to  the  real  axis.  Because  sin  Ax  is  an  odd  function  and  kQ  -X  -  if -v  is 
an  even  function  of  A,  the  integral  is  zero.  The  same  is  true  for  the  pro¬ 
posed  integrations  over  [i  and  v  .  The  integral  then  is: 


U  ~ 


*  C  f  T  CO*  1  a:  «in  2  a  i.  —  «  «in  i  y  wn  2  a  n 

in*  J  J  J  co*  2a  l  —  cos  2  a  a 


dldfid  v 

' '  ’  k*  -  i*  -  fi*  -  »•»  ‘ 


(29) 


It  is  easy  to  convert  this  integral  into  a  sum  built  tip  from  the  residues 
at  the  poles  of  the  integration.  The  poles  for  the  A-integration  are-given 

by 

cos  2  a  ^  «■  cos  2  a  a ,  A0  =  --  ^  a-- , 

where  I  is  a  positive  or  negative  integer.  In  the  A-plane  the  poles  lie  on 
the  two  parallels  to  the  real  axis  which  go  through  the  point  ±a .  (Figure  4). 

The  integrand  is  an  odd  function  of  A.  Then  by  the  substitution 
A*  -  -A,  the  horizontal  original  path  (1)  is  replaced  by  path  (2)  which  lies 
symmetrically  tc  (1)  below  the  lower  row  of  poles  and  is  traversed  from 
■fcoo  to  -oo.  Adding  a  factor  1/2  for  integration  in  each  variable,  we  may 
use  the  path  (1)  +  (2),  thus  obtaining 
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so*  1  x  «in  2  a  l  -  *  *i&  A  x  tin  2  a  a 
cos  i  a  A  —  cos  2  a  a 

didfidr 


“'k*  -  i*  -  fi*  -  *S1 


Tbe  path  (1)  +  (2)  is  a  path  closed  at  infinity  which  can  be  contracted  to 
circulation  of  the  row  of  poles,  the  circulation  occurring  in  the  negative 
sense. 

The  residue  at  one  cf  the  poles  X^  =  (/  x±  ao)/a  is  obtained  by  develop¬ 
ment  of  the  integrand  denominator  in  terms  of  (X  ~  ^q)  : 


co*  Iq  x  c:o  2  g  1«  —  «  tin  ^  a:  sin  2  a  o 
—  2  n  sio  2  a  i,  (1  —  !«}  4-  . . . 
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where  the  upper  signs  refer  to  the  poles  of  the  upper  set  and  conversely. 
Tbe  sum  of  all  residues  which  result  from  the  integration  over  X  is 


or  since  the  second  sum  is  identical  with  the  first,  (differing  only  by 
having  -i  in  place  of  i  for  the  summation  index): 


By  the  same  procedure  regarding  the  other  integrations  we  obtain  the  total 
potential: 
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(30) 


This  series  which  is  equivalent  to  the  grals  Eqs.  (28)  and  (29)  is  a 
precise  description  of  the  electromagnetic  vector  potential  which  results 
from  the  superposition  of  the  spherical  waves. 

By  placing  the  factor 

£  -*(•*  +  /» +  r«)  B  £  —  ik» 

before  the  summation  sign,  the  potential  has  the  form 


77  « 


n 

tube 


(30a) 


where  the  remaining  sum  E  is  a  periodic  function  (Fourier  series)  with 

"“iks 

the  development  range  ±  a,  ±b,  ±c.  The  factor  e  differs  from  the 
factor  for  the  dipole  phase  only  therein,  that  in  place  of  the  discontinuous 
variable  S  (the  coordinate  of  a  dipole)  the  continuous  coordinate  s  is  sub¬ 
stituted. 

In  order  to  obtain  a  clear  picture  of  the  potential  given  by  Eq.  (30a), 

one  must  keep  in  mind,  that  the  development  range  of  the  periodic  function 

is  extremely  small  in  comparison  with  the  distance  that  the  factor  e 

requires  for  a  full  period.  The  distance  between  the  atoms  in  the  body  are 

of  tbe  order  of  10  cm  while  a  light  wave,  even  in  strongly  refractive 

-5 

material,  is  not  smaller  than  10  cm.  There  are  thus  100  -  1000  atoms 
to  a  light  wave  length. 
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The  factor  e-1*  S  thus  hardly  changes  within  an  elementary  paralle- 
piped.  In  order  to  visualize  Eq.  (30a)  we  have  primarily  to  think  of  a 
periodic  function  of  very  complicated  character  repeating  around  each  atom. 
After  that,  proceeding  from  atom  to  atom,  we  have  to  modify  the  ampli¬ 
tude  of  the  function  according  to  the  modulation  given  by  the  slowly  changing 
factor. 

For  the  optical  wave,  roughness  and  spikes  of  the  microscopic  field 
are  immaterial  and  the  periodic  function,  in  its  repetition  from  atom  to 
atom,  contributes  only  by  its  average  value.  This  is  given  by  the  constant 
term  000  of  the  Fourier  series,  namely 


1 

V-** 

Therefore  the  potential  of  the  optical  field  is: 


q.-_ e 

t*be 


*»- V* 


(31) 


completed  by  the  time  factor  e  ^  it  describes  a  plane  wave  in  the  direc¬ 
tion  -s  whose  exponent  contains  k  as  the  wave  constant  and  not  kp  as  do  the 
individual  spherical  waves.  Accordingly,  the  optical  potential  satisfies 
the  wave  equation 


An  +  k*n  =  o. 


The  wave  constant  k,  which  was  originally  assumed,  and  the  ratio 

k 


have  to  be  determined  from  the  "dynamic  equilibrium  condition"  (Eq.  10, 
Section  2) . 
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The  factor  (k0?<  -  k2)  1  dete mines  the  amplitude  a  of  the  oscillating 
electrons  for  a  given  intensity  of  the  observed  optical  wave.  For  small 
values  of  (k^  -  k),  namely  when  the  index  of  refraction  is  not  far  from  one, 
the  amplitudes  are  particularly  small.  This  is  fundamentally  the  reverse 
of  the  statement,  that  slightly  resonating  electrons  have  only  a  small  in¬ 
fluence  on  the  index  of  refraction. 


7.  SECOND  SUMMATION:  EXCITATION  FIELD 

The  form  of  the  total  potential  as  derived  in  the  previous  paragraphs, 
shows  clearly  the  origin  of  the  optical  field.  In  particular,  it  was  easy  to 
see  that  the  total  potential  has  the  same  phase  e  at  all  lattice  points  as 
that  of  the  dipoles.  This  means  that  in  the  calculation  of  dipole  oscillation 
(Eq.  10,  Section  2)  we  may  restrict  ourselves  to  the  dipole  at  the  origin. 
For  it  we  have  to  find  the  excitation  potential,  or  rather  its  second  deriva¬ 
tive,  which  by  means  of  the  determinant  (Eq.  11,  Section  2),  determines 
the  phase  constant  k. 

It  is  not  possible  to  separate  the  effect  of  the  test  dipole  000  from  the 
sum  Eq.  (30)  of  the  previous  paragraph.  This  is  due  to  the  fact  that  the 
sum  becomes  infinite  for  x  =  y  =  z=  0by  non-uniform  convergence.  The 
contribution  which  has  to  be  subtracted,  R  *  e^^  ’  diverges  in  an  entirely 
different  manner.  It  is  not  possible  to  combine  both  terms  before  they 
become  infinite.  Furthermore,  the  form  in  which  II  appears  does  not  allow 
dou^.e  differentiation. 

fi  is  very  fortunate,  especially  with  regard  to  the  numerical  calcula¬ 
tions,  that  the  summation  also  may  be  performed  with  the  help  of  the  single 
integral  for  the  individual  potential  [Eq.  » 18),  Section  5],  One  obtains  a 
single  integral  over  a  product  of  theta  functions  and  exponential  functions 
and  one  can  effectively  perform  the  subtraction  of  the  divergent  terms  from 
the  integrand  so  that  in  this  calculation  of  the  excitation  potential,  no  in¬ 
definite  term  »  -  cc  occurs.  This  is  fundamentally  the  merit  of  the 
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d  -transformation  theorem,  by  virtue  of  which  each  0  -function  admits  of 
two  entirely  different  developments.  In  one  development,  each  term  can  be 
interpreted  by  us  as  the  effect  of  a  dipole,  and  therefore  the  contribution  of 
the  test  dipole  to  the  total  potential  can  be  recognized  and  subtracted.  In 
the  other  development,  the  contributions  cf  the  individual  dipoles  are  com¬ 
pletely  mixed  up  and  each  term  of  the  theta-expansion  matches  one  of  the 
terms  from  which  in  the  previous  sections  we  have  built  up  the  total  poten¬ 
tial  in  the  form  of  the  triple  sum.  By  suitable  use  of  both  descriptions,  it 
is  not  difficult  to  subtract  the  contribution  of  the  test  dipole,  while  retaining 
the  essential  parts  of  the  fields  of  the  other  dipoles. 

We  proceed  from  the  integral  Eq.  (18),  Section  5: 
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in  which  we  substitute  a  dimensionless  variable  e  by  setting 


«.»-**■*  jp»_ 
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where  d  is  a  length  of  the  order  of  magnitude  of  a,  b,  c  which  we  will 
choose  later  as  the  geometric  mean  of  the  lengths  of  these  sides.  Then 
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Hie  condition  of  Section  5,  regarding  the  direction  under  which  the  path  of 
integration  has  to  leave  the  origin  is  not  affected  by  this  transformation. 
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We  multiply  by  tbe  dipole  phase  factor  and  sum:  This  gives  for  the  total 
potential 


n  = 


(Oj 


'j-  (35) 
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The  part  of  the  exponent  which  is  depended  on  the  summation  index  l  is: 


-*r^P  +  (^-2iajl. 


The  simplest  theta  function  (which  is  referred  to  also  as  &m  or  d,  ) 
admits  of  the  following  series  expansion: 


+  • 


(36) 


By  interchanging  the  sequence  cf  summation  and  integration,  we  have  there¬ 
for  in  Eq.  (35)  before  us 

„  /bo  . ax  «j  . a* 

and  the  total  potential  is 
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2  we  desire  to  drop  the  dipole  000  from  the  summation,  we  omit  the  000 
term  and  we  obtain  the  excitation  potential : 


IT  =-  \  ft  &  &  &  -  i;« 


^  dt . 


in  order  to  perform  the  integration,  we  utilize  the  series  expansion  oi 
&  .  The  original  series  £q.  (36)  converges  quickly  when  t  is  large;  we 
obtain  a  useful  expansion  for  small  e  through  tbs  v  -transformation 
formula^' 


&{v  ‘  it) 


(39) 


By  means  of  the  ’transformer?’  and  the  ’brigmal”  series,  the  following  in¬ 
tegrals  for  II  are  produced: 


JL.-r_  J'ST*  ‘  **  1 '  •  *>  *  *  t~3d* 

d  abcj 


The  integrands  are  identical!/  equal  for  each  value  of  e .  We  will  abbreviate 
the  exponents  in  the  first  and  second  integral  as  I  and  11  respectively.  Then 
the  derivatives  are: 


(1)  See  for  example,  A.  Eraser,  Lehrbuch  der  Theta  fur  kticmen, 
p.  96,  Teubner. 
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We  must  now  subtract  from  II  the  second  derivative  with  respect  to  x  of 
the  individual  potential  of  the  test  dipole  at  the  origin*  For  this  purpose, 
we  divide  the  integration  for  II  in  two  parts:  from  0  to  E  and  from  E 
to  x  .  The  crossover  place  E  remains  at  present  indeterminate.  On  the 
first  part  of  the  path  we  use  the  series  with  the  exponent  L  on  the  second, 
the  series  with  IL  The  quantity  to  be  subtracted  is  [  see  Eq.  (34),  this 
Section  ] : 


%tr 


*  d, . 


(43) 


Again,  we  divide  this  integral  in  two  parts  at  the  crossover  point  E  ,  and 
find  that  the  integral  from  E  ro  so  is  identical  with  the  term  000  of  the 
second  half  of  0^.  The  subtraction  can  therefore  be  achieved  by  putting 
an  accent  oo  the  sum  sign  in  order  to  indicate  that  the  term  000  has  to  be 
omitted.  If  we  simultaneously  take  the  value  atx  =  y  =  z  =  P=  0we  obtain 


*  As  published,  due  to  typographical  error,  the  last  power  within  the 
brackets  was  given  as  S. 
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where  the  exponents  at  the  origin  are: 


I 

II 


in  n  ^}  the  general  sun  term  may  easily  be  integrated  when  one  takes  into 
consideration  that  e  "3d£  ^fdd/t2)  .  The  indefinite  integral  produces  no 
contribution  at  f  =  0  because  cf  the  direction  of  the  departure  from  the 
zero  point;  there  remains  the  value  at  the  upper  limit  E  and  one  has: 


Here-  one  sees  plainly  the  relation  to  the  total  potential  in  the  form  of 
Section  6.  For  E  =*,  IT^  goes  over  into  the  (formally)  differentiated 
total  potential  Eq.  (30),  Section  6,  taken  at  the  origin  and  diminished  by  the 
(likewise  infinite)  test  dipole  term,  5  £  is,  however,  given  a  finite  value, 
then  all  the  terms  are  finite  because  the  infinite  effect  of  the  test  dipole 
has  been  taken  out  ex  a*2*. 

XX 
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Through  a  similar  calculation  the  following  is  obtained: 
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In  contrast  to  II  there  exists  no  term  in  II  ^  which  corresponds  to  the 
xx  xy 

single  potential  to  be  subtracted.  The  reason  for  this  is  that  in  the  dif¬ 
ferentiation  of  Eq.  (34)  with  respect  to  x  and  y  no  term  appears  which  be¬ 
comes  infinite  independently  of  the  direction.  In  spite  of  the  factor  i  •  m 

(21 

the  accent  on  the  sum  in  II  ^  may  not  be  omitted,  because  the  integral 
over  the  term  !  =  m  =  n  =  0  diverges. 

The  series  Eq.  (45)  for  11^  and  (44)  for  11^  give  the  necessary  data 

XX  *  ^  'r 
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for  the  calculation  of  the  principal  indices  of  refraction  inasmuch  as  it  is 
possible  to  find  them  without  the  use  of  approximations.  In  the  following 
sections,  we  shall  make  use  of  the  simplifications  resulting  from  the  actual 
orders  of  magnitude  in  the  case  of  propagation  of  visible  light  in  solid 
bodies. 


8.  APPROXIMATIONS  DL!E  TO  THE  SMALLNESS  OF  *,  v 

Equations  (42)  and  (43)  established  in  Section  7  for  the  derivatives  cf 
the  excitation  potential  at  the  place  of  the  dipole  000  can  be  considerably 
simplified  when  x»,x  and  d.ce,  d./3,  d.7  are  small  numbers.  This  is 
always  the  case  when  visible  light  in  solid  material  is  treated,  due  to  the 
fact  that  d  is  of  the  same  order  as  the  interatomic  distances  and 
kg,  k,  o,  y  are  of  the  order  of  l/X,  As  already  mentioned,  the  ratio 
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of  dA  is  of  the  order  of  1/7.00  to  1/1000  and  the  above  quantities  are  of  the 
same  order. 

Under  these  c  rcumstances,  terms  suer  as  da/*  occurring  in  the  series 
Eqs.  (44)  and  (ao)  Section  7,  play  a  role  only  v/hen  i  -  m  =  n  =  0. 

We  omit,  therefore,  in  all  terms  of  the  series,  the  small  values  after 
having  segregated  the  term  000.  This  term  is  in  according  to  Eq.  (14), 
Section  7: 

n*  <Pn*  1  n  co*%  (8,z) 

2abe  n*  *,■'-«*  ""  2 abe  1-V  ’ 

where  *  *=  :-/x0  is  the  index  of  refraction.  The  term  is  therefore  of  the 
order  of  magnitude  1  in  comparison  to  x*  . 

Because  of  the  smaliuess  of  the  integral  appearing  in  Eq.  (45)  of 
the  prec~  section 


je  '  e*de 

2 

has  an  integrand  which  differs  appreciably  from  f  only  where  the  integrand 
is  very  small  anyhow.  We  commit  an  error  at  most  of  the  order  *oa  , 
when  we  replace  this  integral  with  one  indepe  dent  of  Kq  : 


s 


m 


A  lter  these  steps,  the  first  part  of  the  excitation  potential  becomes 
nil  -  TOT  (t r,r  «“■(«,*)  +  *g). 
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where  in  this  approximation  ip\,  depends  neither  on  x  and  nor  on  the 

AA 

direction  of  propagation  of  the  wave,  but  only  on  a,  b,  and  c. 

/  2> 

Similar  considerations  are  applied  to  IP  7  Eq.  (44),  Section  i.  Also 
here  the  small  quantities  of  the  order  of  Xq*  can  be  disregarded  in  the 
exponent  of  the  e -function  wherever  they  stand  next  to  terms  of  the  order 
of  magnitude  1.  Only  for  the  term  000  would  this  not  be  the  case;  but  this 
term  is  excluded  because  of  the  accent.  Thus  II  ^  becomes  entirely  in¬ 
dependent  of  the  small  quantities.  We  write  in  this  approximation 


/7(2)  -_J!L_  ti/(2> 

**  iabe  •  “ 


(1)  (2) 

and  obtain  altogether  for  the  excitation  potential,  after  setting  ip  +  \p  =  </ 


H'.z  “  ~2abe'  (t=7T-  C0S*  ^ 


and  similarly 


-  w(t^-  —  « *>  “» «  jo)  • 


In  n’y  and  analogous  derivatives,  it  is  evident  that  no  term  ip^  appears; 
for  this  term  would  be  independent  of  the  direction  of  propagation  and  thus 
would  also  exist  in  che  case  of  propagation  along  one  of  the  coordinate 
axes.  But  we  know  from  symmetry  that  in  this  case  IT^  =  0  (see  above) . 
Finally,  we  have  in  this  notation: 
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9.  THE  DISPERSION  FORMULA 

In  Section  2  we  had  to  distinguish  six  cases  of  principal  indices  of  re¬ 
fraction,  depending  on  the  directions  of  the  oscillation  and  of  the  ray.  But 
we  have  previously  suspected  that  pairs  of  these  indices  would  so  nearly 
coincide,  that  their  discrimination  would  be  practically  superfluous.  As  a 
matter  of  fact,  we  can  now  write  down  the  equation  for  the  propagation  of 
the  optical  wave  in  the  z  direction  with  oscillation  in  the  x  direction: 


fl  +(JZ7'-  /TJ..J.J  -  12  +  ^(j-^  +  ~  V..)  -  0 


and  we  find,  that  in  this  approximation,  because  of  the  independence  of  ip 

from  the  direction  of  the  ray,  the  identical  equation  would  be  valid  for  v  , 

that  is  for  a  =  7=0,  /3  =  k  .  We  therefore  omit  the  distinction  between 

and  v _ ,  and  we  write  v„  for  the  principal  index  of  refraction  for  dipole 

XV  •  1  -  ■  ■  xz  .  X - - 

oscillations  parallel  to  the  x  axis. 

Let  v  be  the  volume  8abc  which  contains  one  dipole;  then  for  the  index 
of  refraction  we  have  the  equation 


*,*-  1 


—-Q  +  A  ip  —  \pa 


(50) 


or 


1  -f 


—  Si  +  a  y  -  - : 


(51) 


We  recall  the  significance  of  the  electron  constant  f2  [Eq.  (9), 
Section  2] 


Si 


-3Hv  s); 


For  a  transparent  crystal  we  set  the  frictional  coefficient  g  =  0 .  Then 
the  dispersion  formula  is: 


v*  =  1  + 


1  + 


l 


t  m  .  .  .  .  v  m  t 

-- V-(i + 


B,  —  An 


*  j 


(52) 

(53) 


or,  introducing  the  wavelength  X  = 


2£c 

n 


=  - - - 

•  T  1*  A,  -  4  n1 


C*A 


(54) 


We  have  here  before  us  a  dispersion  formula  with  tv/o  constants,  of  the 
sane  type  as  those  of  H.  A.  lorentz*  and  of  Planck.  The  latter  would  be 
in  our  terminology 


v3  =  1  + 


v  m  t  1  cm, 
Im*"0  3  4 lie' " 


whereas  Lorentz  has  -(1/3  +  s)  in  place  of  - 1/3  where  s  remains  indeter¬ 
minate.  In  all  three  formulas  this  term  indicates  the  influence  of  the  di¬ 
poles  of  the  body  on  a  test  dipole.  The  terra  1/3  arises  from  the  contri¬ 
butions  of  the  more  distant  surroundings  of  the  test  dipole,  where  the  di¬ 
poles  may  be  spread  out  into  a  homogeneous  medium.  In  Lorentz'  tern  s 
is  added  the  strong  but  unpredictable  influence  that  near  neighbor  dipoles 
exert  on  the  test  dipole.  For  a  cubic  arrangement  of  neighboring  molecules, 
H.  A.  Lorentz  proves  that  s  =  0,  so  that  in  this  case  his  formula  becomes 
identical  with  Planck's. 


*Cited  under  Section  4 
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Also  from  our  equations  specialized  for  a  cubic  lattice  the  term  1/3 
is  obtained.  For  in  that  case  one  has,  by  symmetry 

1  —  Axft  +  tft,,  ■«  1  —  |  Ay.  (55) 

But  A \i/  can  be  determined  in  the  following  maimer:  one  develops  the  total 
potential  Eq.  (30),  Section  6  in  powers  of  k  obtaining 

when  only  the  lowest  powers  are  written  out. 

In  order  to  obtain  the  excitation  potential. 


i 


should  be  subtracted  from  this.  The  term  with  k  in  the  development  of 
II  goes  over  unchanged  in  the  development  of  O'.  Thus  one  has 

ah- -t,‘n- +  +  •••]) 

-v(i=?-  +  C*-  +  -J) 

“v(i^  +  I  +  £*•••])• 

By  comparison  with  Eq.  (49)  Section  8  it  follows  that 


Ay  —  1. 

and  thus  from  Eq.  (55) 


1  -  Ay  +  ym  = 
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Insofar  as  the  isotropic  body  may  be  replaced  by  a  cubic  arrangement, 
our  formula  thus  entirely  agrees  with  the  Lorentz-PIanck  formula. 

Discussion  of  the  dispersion  formula  thus  obtained  would  add  nothing 
new  regarding  the  dependence  on  the  wavelength  of  tne  indices  of  refrac¬ 
tion  and  absorption  than  what  has  already  been  discussed  by  Planck. 

Only  those  details  which  are  produced  through  the  anisotropy  of 
Eq.  (52)  shall  be  noted  here: 

(1)  From  two  measurements  of  v  .  the  quantities  A  and  B  [  Eq.  (53)  ] 
can  be  computed.  The  quantities  o  are  functions  only  of  the  lattice  constants 
a,  b,  c.  Thus  by  using  the  value  of  a  single  principal  index  of  refraction 

at  two  different  wave  lengths  one  can  obtain  the  complete  dispersion  formula 
not  only  for  this  index,  but  also  for  the  two  other  principal  indices  of  re¬ 
fraction. 

(2)  The  constant  A,  independent  of  which  index  of  refraction  would  be 

2 

used  in  its  determination  should  have  the  value  i//4'*m/e  which  is 
characteristic  for  the  lattice.  The  same  value  of  A  should  therefore  be 
found,  no  matter  fiom  which  principal  index  of  refraction  it  is  obtained. 

But  one  finds  this  relationship  not  fulfilled.  For  example  for  anhydrite, 
(anhydrous  calcium  sulphate) 

=  1.6125- 10-33;  Ap  =  1.6726-10"”;  Ar  =  1.5501  -lO"3', 

values  which  deviate  by  8  percent. 

One  of  the  reasons  for  this  deviation  is  the  fact  that  the  electrons 
having  the  proper  frequency  n^  do  not  form  the  simple  Bravais  lattice  that 
has  been  presupposed  here.  A  second  reason  is  that  in  addition  to  the 
electrons  with  the  characteristic  frequency  nQ  others  with  different  eigen 
frequencies  and  possibly  ions  with  eigen  frequencies  in  the  infrared  are 
contained  in  the  crystal. 

The  two-constant  Eq.  (52)  is  then  to  be  understood  as  an  approximation 
to  a  formula  of  complicated  structure  and  the  constant  A  of  the  dispersion 
formula  conse  quently  loses  its  simple  meaning  and  can  vary  with  the 


! 
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direction  of  the  propagation  of  the  wave  or  of  the  polarization.  In  this 
connection  the  observation  of  Drude^  is  relevant:  if  in  calcite  incidence 
is  at  right  angles  to  the  optical  axis,  the  ordinary  ray  involved  apparently 
five  electrons,  and  the  extraordinary’  ray  two  electrons  per  molecule. 

(3)  From  Eq.  (52)  it  follows  immediately  that  the  effect  of  the  lattice 
arrangement  can  be  interpreted  also  as  an  anisotropic  bonding  of  the 
electrons  in  the  molecules.  The  variation  in  the  formula  of  dispersion  in 
passing  from  i  ^  to  or  from  to  p  ,  respectively  can  be  replaced, 
as  Eq.  (52)  shows,  for  all  wavelengths  bv  a  suitable  modification  of  the 
eigen  frequency  nQ  .  Located  within  the  lattice,  a  dipole,  whien  in  itseli  is 
bound  isotropically,  acquires  through  the  influence  of  the  structure  three 
different  eigen  frequencies.  The  two  views,  namely  attributing  the  origin 
of  birefringence  either  to  the  inherent  anisotropy  of  the  ultimate  particles, 
or  to  the  anisotropic  arrangement  of  isotropic  particles,  thus  cannot  differ 
in  their  qualitative  consequences. 

10.  THE  EQUATION  OF  THF  NORMAL  SURFACE 

We  have  dealt  up  to  now  exclusively  with  the  principal  indices  of  re¬ 
fraction,  six  of  which  we  had  to  distinguish  originally.  With  neglect  of  the 
terms  of  the  order  (a/A)  we  recognized  that  indices  referring  to  the  same 
polarization  are  equal  to  each  other.  This  made  the  properties  of  the  lat¬ 
tice  equal  to  those  of  a  bi-axxal  crystal. 

The  question  is  whether  this  equivalence  is  true  for  all  directions,  that 
is,  whether  the  relationship  between  the  direction  of  the  wave  normal  S  and 
the  index  of  refraction,  which  is  contained  in  the  determinant  (Eq.  11, 
Section  2),  is  the  same  as  in  the  orthorhombic  crystal. 

Retaining  the  same  degree  of  accuracy,  we  introduce  the  values  of 
Eqs.  (47),  (48),  and  (49)  from  Section  8  in  the  first  line  of  the  determinant. 

Ann,  d.  Phys.  14:691,  1904. 


/I' 
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This  line  becomes 

"  "**  tabc  (f  ITT  ( 1  —  cosI  (*» z))  +  drp—  u»_)  ; 

n  _ _  v  J.  J.; 

2oie  1-**  003  cos  s>y)-  ~  Yabc  °°S  cos  -*c* 
We  introduce  the  abbreviations: 
jV  =  -j~l  <  cos  (*>  *)  -  cos  1  . . . 

and 


2ai>e  ^  ,  ,  >  i 

T  -  £  +  A  xi>  -  xpxx  ==  =  a; 


Thus  the  determinant  reads: 


-V  sin*  1  -  Nx 

—  N  cosl  cos  2 

—  N  cos  1  cos  3 


-  A*  cos  1  cos  2  -  A*  cos  1  cos  8 
A7 sin*  2  -  X,  -  A*  cos  2  cos  8 

—  A'  cos  2  cos  8  AT  sin2  8  —  Nx 


=  0. 


We  write  in  the  diagonal  (1  -  cos2. . . )  in  place  of  sin2. . 
the  factor  which  is  non- zero  except  in  trivial  cases: 

—  N 3  cos2 1  •  cos*  2  •  cos2  8 

There  remains 


.  and  take  out 


(i- A-dM 

{  X  cos’  1  / 


\  iVco*’  2 


(l-Xz*.) 

\  x  COB*  3/ 


=  0 


ila_ 
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or,  except  for  a  factor: 

IF  cot*  1  ,  X  cot'  2  X  cot*  S  .  n 

-jr: xr  +  ~y~~xT  +  ~x^x. —  u* 

This  equation  is  nothing  other  than  that  of  the  normal  surface  of  a  bi-ax ial 

o  2  2 

crystal.  Replacing  the  1  by  cos“l  +  cos“2  cos  3  and  the  N. . .  by  their 
egression  in  v,  one  obtains  the  equation  in  its  usual  form: 

co*’  1  ,  cos*  2  ,  r,*  ec*’  3  _  n 

+  +  >*-*,'  ~ 

Thus  we  see,  that  a  point  lattice  of  the  kind  considered,  behaves  with 
respect  to  all  directions  as  a  crystal  and  that  from  the  s’andooint  of  optics 
there  is  no  objection  to  even  the  simplest  lattice  mode-  of  a  crystal. 


li.  THE  BIREFRINGENCE  D 

Already  the  general  considerations  of  Section  2  have  indicated  a 
relationship  between  the  principal  indices  of  refraction,  wnich  according 
to  Eq.  (14)  Section  2,  is 


{A  IT—  n'xz\  =  ( Air-ir ^ .  (56) 

The  indices  z  on  the  parentheses  indicate  that  on  both  sides  of  the 
equation  the  potentials  are  those  for  the  direction  of  propagation  along  z, 
but  that  on  the  lefthand  side,  prlarization  is  along  x  and  the  refractive 
index  is  ux,  while  on  the  righthand  side  polarization  is  along  y  and  the 
refractive  index  is  v  . 

If  we  insert  in  this  relationship  the  approximate  values  of  Fqs.  (47), 
(48),  and  (49)  from  Section  8,  there  results 


<s  be 


or  since  At  is  the  same  in  both  cases,  we  obtain 


0 


-rA~’jl — «  *  p  -JD  . 

-r/Xl-O  •"  '  »»  r» 


(57) 


The  dimensionless  quantity  on  the  leit,  which  was  recognized  already  by 
T.  H.  Havelock  as  independent  of  the  wavelength,  we  find  dependent  only  on 
the  ratios  a:b:c  of  the  lattice  periods.  We  call  it  the  measure  of  structural 
birefringence. 

Through  the  series  Section  7,  we  are  in  position  to  compute  these 
constants  for  arbitrary  ratios  a:b:c.  For  a  suitably  chosen  crossover 
point  E ,  the  series  converge  rapidly  enough  to  make  the  numerical  calcula¬ 
te  ->n  easy. 

But  at  present  it  is  not  yet  premising  to  undertake  the  calculation  for  a 
known  crystal  in  the  hope  of  finding  agreement  between  calculated  and  ob¬ 
served  values.  None  of  the  crystals  whose  structure  is  known  up  to  now 
(Fall  1915)  has  a  structure  so  simple  that  one  could  assume  with  some 
measure  of  certainty  that  a  single  kind  of  electron  is  adequate  to  render 
its  optical  properties. 

The  numeric?  1  computation  and  the  extension  of  the  theory  to  more  com¬ 
plicated  spsiee  lattices  shall  therefore  be  omitted  here. 

Only  in  order  to  give  a  general  idea  about  the  expected  order  of  the 
magnitude  of  the  effects,  we  take  for  the  calculation  the  ratios 
a:b:c  =  0. 8932  :  1  :  1. 008  which  are  the  same  as  those  of  anhydrite.  Cal¬ 
culation  leads  to 

DXf  -  0.17848  JOxt  =  0.17983  Dyt  =  0.00133, 
whereas  the  observed  values  for  the  same  material  are: 


D}.  =  0.05145  I>Ta 


G.05925  JDy/  =  0.00854 


From  this  it  may  be  concluded  that  the  anisotropic  arrangement  alone  is 
sufficient  for  explaining  the  order  of  magnitude  of  the  observed  birefringence. 

That  the  simple  assumptions  concerning  the  space  lattice  do  not  suffice 
for  rendering  the  phenomena,  is  also  indicated  by  the  fact  that  the  observed 
values  of  D  are  found  independent  of  the  wavelength  only  to  a  moderate 
extent.  Thus  in  quartz  the  variation  of  JD  between  the  lines  C  and  E  of  the 
spectrum  is  15  times  smaller  than  the  variation  of  tbs  double  refraction 
(u  -e)  in  the  customary  sense.  In  most  crystals  the  variation  of  D  is 
some  3  to  6  times  smaller  than  that  of  the  usual  birefringence,  but  there  are 
crystals  for  which  D  varies  more  than  the  birefringence. 

Whether  the  constancy  of  D  could  be  valid  for  arbitrarily  complex  point 
systems  constructed  of  one  kind  of  electrons  may  remain  undecided.  Cer¬ 
tainly  it  is  not  valid  for  complex  point  systems  with  different  sorts  of  ions, 
since  it  is  not  possible  to  eliminate  the  two  or  more  electron  constants  n. 
from  the  two  dispersion  equations  for  v  and  v  . 

^  y 

Unfortunately  no  similar  simple  relationship  for  two  sorts  of  electrons 
or  ions  could  be  found  between  the  three  principal  indices  of  refraction, 
which  would  be  valid  regardless  of  wavelength. 
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Part  II.  Theory  of  Reflection  and  Refraction* 


1.  INTRODUCTION 
LI  Problem 

Let  the  upper  half  space  be  filled  with  a  lattice  of  dipoles  with  the  co¬ 
ordinates 

.Y  =  2 ol;  Y  --  2 bm;  Z  —  2cn.  (58) 

j 

! 

where  i,  m,  n  are  integers  and  i,  m  assume  all  positive  and  negative  values, 
while  n  assumes  only  positive  values  (including  0).  Further  let  a  plane 
wave  be  incident  on  the  lattice  from  below  and  let  its  wave  normal  have  the 
direction  (see  Figure  1).  Its  periodic  function  is 


Find  the  mode  of  oscillation  of  the  lattice  and  the  thereby  created 
optical  field. 


*  Continuation  of  the  work  by  the  same  title  in  Ann.  d.  Phys.  49:1-38, 
1916,  quoted  as  Part  L  "" 


interior 


Exterior 


Figure  5 


1.2  Procedure  of  Computation 

A  procedure  which  appears  to  hold  little  promise  would  be  to  write 
down  methodically  the  excitation  equations  for  all  dipoles  and  to  solve  them 
beginning  with  those  at  the  boundary.  The  force  terms  appearing  on  the 
right  hand  side  of  these  equations  would  themselves  depend  on  all  of  the 
dipole  oscillations  and  at  that,  in  a  most  complicated  manner.  Consequently 
the  computation  of  the  force  without  previous  knowledge  of  the  dipole  oscil¬ 
lations  is  by  no  means  practical. 

Another  method  of  solution  exists  in  making  a  suitable  assumption  for 
the  dipole  oscillations.  We  can  then  actually  sum  the  electric  forces  and 
satisfy  ourselves  that  the  dynamic  conditions  are  fulfilled  for  the  assumed 
mode  of  oscillations. 

Therefore,  the  assumption  is  made  here  that  in  the  crystal,  up  to  the 
very  boundary,  there  exist  planes  of  equal  dipole  phase.  In  the  unbounded 
crystal  such  an  assumption  leads  to  the  propagation  of  a  plane  optical  wave. 
Let  the  wave  normal  of  this  dipole  wave  be  s,  and  its  velocity  q. 

Then  the  summation  of  the  spherical  waves  issuing  from  the  dipoles 
does  not  yield,  as  in  the  case  of  the  unbounded  crystal,  only  the  field  of  a 
plane  wave,  which  propagates  in  the  same  manner  (velocity  q,  direction  s>/ 
as  the  dipole  phase.  Rather,  there  occurs  besides  this  field  inside  the 
crystal,  a  plane  wave  with  the  velocity  c  whose  direction  of  propagation  sfj 
is  related  to  s  by  the  law  of  refraction.  Moreover,  in  outside  space,  a 
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wave  comes  into  existence  wMch  is  a  mirror  image  with  regard  to  the  sur¬ 
face  of  the  last-mentioned  internal  wave  and  it  has  likewise  the  velocity  c. 
The  outside  wave  is  the  reflected  wave.  The  first  inside  wave  (q)  is  the 
refracted  wave.  The  wave  on  the  inside  of  velocity  c  is  the  one  which  is 
required  by  dispersion  theory  in  order  to  cancel  the  incident  wave.  It 
begins  with  full  strength  at  the  boundary  and  can  therefore  compensate  the 
incident  wave  Sq  throughout  the  interior  of  the  crystal.  (It  is  called  for 
short,  the  internal  boundary  wave. ) 

This  result  already  offers  the  solution  of  the  dynamical  conditions  for 
the  dipole  oscillations.  For  when  this  same  boundary  wave  and  the  incident 
wa?e  cancel  each  other  inside  the  crystal,  then  only  the  same  field  remains 
rnside  the  bounded  crystal  as  in  the  unbounded  crystal,  and  we  only  have  to 
see  to  it  that  the  assumed  dipole  oscillations  are  compatible  with  this  field. 
This  is  accomplished  when  the  propagation  velocity  of  the  refracted  wave 
has  the  correct  value,  as  was  shown  in  the  dispersion  theory. 

Amplitude  and  phase  of  the  internal  boundary  wave  are  determined  by 
the  requirement  that  the  incident  wave  be  cancelled.  Since  the  three  waves 
which  are  generated  by  the  dipole  oscillations  show  definite  amplitude  and 
phase  relationships,  the  amplitudes  and  phases  of  the  reflected  and  re¬ 
fracted  waves  are  also  determined  by  this  condition  of  cancellation.  The 
resulting  relationship  between  the  amplitudes  of  the  three  waves  must  be 
that  given  by  the  Fresnel  formulas  provided  that  the  correct  transition  is 
made  from  the  microscopic  field  to  the  phenomenological  Maxwell  field. 

Lei  it  be  said  at  once  that  we  restrict  our  thinking  to  the  case  of  a 
cubic  crystal  which  is  thus  also  opt  caiiy  isotropic.  In  this  way  complica¬ 
tions  will  be  avoided,  which  wou'd  only  detract  from  the  essentials  of  the 
theory.  Notwithstanding  this,  we  distinguish  in  the  formulas,  .the  three 
axis  intercepts  a,  b,  c,  which,  if  desired,  immediately  may  be  considered 
as  equal. 


2.  THE  SUMMATION  OF  THE  DIPOLE  WAVES 


As  in  the  first  part.  Dispersion  Theory,  we  employ  the  Hertz  potential 
of  a  dipole  in  the  form  of  expression  Eq.  (23),  Part  I,  Section  5,  namely 
with  omission  of  the  factor  eo*  e 


£«'*»/•  1  Crr  »•  ±  < ■)  dl  dud ¥ 
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(00) 


(X,  Y,  Z)  are  the  coordinates  of  the  dipole,  (x,  y,  z)  those  of  the  field  point. 
The  variables  A,^,  v  run  parallel  to  the  real  axis  between  it  and  the  point 
kQ,  (which,  initially,  is  thought  of  as  having  a  positive  imaginary  part) 
(Figure  6). 

Tfco  potentials  to  be  summed  are: 


f 


(61, 


Figure  6 
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Here  (in  the  terminology  of  Part  I,  Section  10) 

k  S  —  /<•  cosl  •  A"  -f  A'  c  -  *  V  4-  k  cos8  •  Z 

=  aX  .  #Y  T  (62) 

This  corresponds  to  a  progressive  excitation  along  +s  . 

On  performing  the  summation  under  the  integral  sign  in  Eq.  (60),  the 
threefold  sum  is  decomposed  into  the  product  of  three  simple  sums.  Of 
these,  the  first  two  whose  indices  l  and  m  run  from  -» to  +»,  are  identical 
with  the  sums  over  /  and  m  which  appear  in  the  corresponding  problem  of 
the  unbounded  crystal  (see  Part  I,  Section  6) .  Our  present  assumption 
differs  only  insignificantly  from  that  made  there,  namely  in  the  signs  of 
O',  p,y.  These  are  easily  changed  in  the  results  of  the  summations.  The 
sum  over  i  thus  produces  here  in  the  integrand  the  factor 
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and  similarly  for  the  sum  over  m. 

On  the  other  hand,  the  summation  over  the  index  n  extends  only  from 
0  to  oo  and  we  have: 


l ’ 


(64) 


hi  order  that  the  sum  converge,  v  +  y  must  have  a  positive  imaginary 
part;  that  is,  v  must  be  conducted  on  a  path  above  -  y  .  In  the  former 
sums,  as  also  in  the  unbounded  crystal,  it  was  required  that  X  be  led 
above  a,  \i  above  3 .  The  feasibility  of  this  requirement  provided  the  sums 
are  unconditionally  convergent  has  been  discussed  in  the  thesis  and  has 
been  attributed  to  the  fact  that  ICq  has  a  larger  imaginary  part  than  has  k. 
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This  assumption  is  also  here  necessary  for  the  summation.  In  the  result 
we  go  without  further  discussion  to  real  k,  (see  Part  I,  Section  3). 

The  total  potential  of  the  "half  crystal"  consequently  becomes: 
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From  this  description,  which  is  valid  for  all  v:-  .u-jS  of  (x,y,  z)  we  pro¬ 
ceed  to  a  description  of  II  as  a  sum  by  reducing  the  Ihreo  integrations  to 
circulations  of  the  poles  of  the  integrand  in  the  ccmple-:  plane  and  evaluating 
these  by  residues.  This  entails  oi  itself  a  distinction  of  inner  and  outer 
space  with  each  containing  different  groups  of  waves. 

The  first  two  integrations  over  X  and  n  are  identical  with  those  in  the 
expression  Eq.  (23),  Part  I,  Section  6  for  the  potential  of  the  v  bounded 
crystal.  They  will  be  treated  just  as  there.  By  utilizing  the  fact  that  the 
integrand  is  an  odd  function  of  the  variable,  the  path  of  Figure  6  can  be 
completed  by  one  naming  below  the  real  axis.  The  series  of  polos  of  the 
integrand  now  lies  between  the  two  paths  at  the  positions 


Residue  development  in  the  integral  over  i  yields  the  sum 
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and  the  same  treatment  is  applied  to  the  integration  over  p.  (Compare 
Parti,  Section  6.) 

By  tin's  transformation  of  the  integrations  over  X  and  p  the  potential  of 
the  bounded  crystal  becomes: 


■T  OO 


This  is  a  mixed  form  of  the  potential,  inasmuch  as  two  integral  signs  have 
been  replaced  by  summation  signs,  while  the  third  integration  still  re¬ 
mains.  This  form  shares  with  Eq.  (65)  the  advantage  that  it  is  valid  for 
the  entire  space. 

We  now  proceed  to  replace  the  last  integral  also  by  a  sum  of  residues. 

To  this  end  we  observe  that  the  integrand  has  two  classes  of  poks  ac¬ 
cording  to  whether  the  first  or  second  denominator  is  zero. 

<0  1  — <.* **;<>  + r>  —  0;  r  =  —  y,  where  n  is  an  integer  (65) 

C  _ 

This  series  of  poles  corresponds  completely  io  the  poles  Eq.  (6^  in  the 
X  and  p  plane,  at  which  the  residues  were  taken.  The  poles  lie  ca  a  straight 
line  parallel  to  the  real  axis  and  below  t,  out  above  the  point  -kp.  This 
line  thus  is  between  the  integration  path  for  v  and  the  point  -k^  (Figure  7). 
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*  The  position  of  this  line  was  stated  incorrectly  when  published  and 
is  corrected  here. 

**  The  words  'that  is’1  have  been  supplied  in  translation  in  order  to 
clearly  separate  the  two  statements  of  Eq.  (70). 


Figure  7 


By  -f  l/~  we  understand  the  value  of  the  root  which  is  situated  in  the  first 

quadrant.  These  poles  form  a  chain  which  sheathes  the  branch  cut  of  the 

* 

root  extending  from  to  infinity  .  With  increasing  i  and  m  the  poles 
become  ever  more  imaginary. 

(c)  Also  the  other  root  value,  thus  under  retention  of  the  positive 
sign  for  the  principal  value,  the  value 

V  =  —  rim 

causes  the  denominator  to  vanish. 

Altogether:  three  lines  of  poles  exist  in  the  v  plane  as  shown  in 
Figure  7. 

The  integration  path  for  v  runs,  as  mentioned,  parallel  to  the  real 
axis  between  the  pole  series  (a)  and  the  series  (b).  This  path  can  either 
be  deformed  upwards  and  then  contracted  about  the  pole  series  (b)  with 
circulation  about  each  pole  in  the  positive  direction:  or  it  can  be  polled 
downward,  namely  in  case  that  the  integrand  vanishes  exponentially  at 

*  Clearly  this  branch  cut  does  not  exist  for  the  integration  over  v 
since  in  the  denominator  only  appears. 
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infinity  in  the  lower  half  plane.  In  doing  so  it  remains  held  by  the  pole 
series  (a)  and  (c),  (both  with  negative  circulation  of  the  poles). 

The  condition  for  the  upward  deformation  is  that  the  factor 


i  _  «-*><»  +  r< 

vanish  at  infinity  in  the  upper  half  plane,  thus 
z  <  0 

This  is  fulfilled  in  all  outer  space. 

The  downward  deformation  on  the  other  hand  is  permissible  because 
of  the  same  factor,  provided  the  denominator  becomes  infinite  more 
rapidly  than  does  the  numerator.  That  is  when 

z  -f-  2c  >  0 

which  is  true  in  the  entire  inside  space. 

These  limits  show  that  a  domain  exists  at  the  boundary 

—  2c  <  z  <  0 

in  which  both  deformations  are  permitted.  The  boundary  domain  continues 
the  crystal  beyond  the  cut  ofi  at  z  =  0  into  outer  space  and  establishes  thus 
the  transition  from  inner  to  outer  space. 

We  have  now  to  obtain  the  residues. 


2.1  Outside  Spacr  -  Negative:  Pole  Senes  so* 


We  have 
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and  thus  the  sum  of  the  residues: 


-idn  —aa)z  a  — 


_ 1 _ 
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(71) 


This  is  the  potential  II  „  which  is  valid  in  the  outer  space.  It  is  a  double 

a 

sum,  wherein  each  term  appears  to  have  the  form  of  an  ordinary  plane 
wave.  In  fact,  this  is  true  only  for  the  term  00  given  that  for  visible  light 
kg,  a,  3,7  are  small  compared  with  l/a,  1/b,  1/c.  For  the  same  reason 
a1!  other  terms  for  which  is  positive  imaginary  represent  inhomo¬ 
geneous  waves.  (See,  however.  Section  7). 

2.2  Inside  Space.  *  4-  2c  Positi\e:  Pole  Senes  (a)  And  (c) 

The  residue  at  a  pole  of  the  series  (a)  is  found  by  development  of  the 
integrand  in  proximity  to  the  pole.  There 

»j.i— cr 

r=ra=:____ 
and  nearby 
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Because  of  the  negative  sense  of  circulation,  the  residue  is  multiplied  by 
the  factor  -2zi.  The  potential  becomes  —  we  designate  it  0^: 
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Except  for  the  already  frequently  mentioned  change  of  the  signs  of  a,  p,  7 
this  is  identical  with  the  potential  of  the  unbounded  crystal.  [  See  Part  I, 
Section  6,  Eq.  (30)  ] . 

Finally,  as  a  second  part  of  the  inner  potential,  the  double  sum  over 
the  pole  series  (c)  is  to  be  added;  but  for  the  sign  of  and  a  further  (-i) 
because  of  the  negative  circulation  this  is  equal  to  the  sum  of  the  residues 
at  the  poles  (b) : 


t;i 
2  ah 


n.= 

22 


(72) 


What  has  been  said  about  the  character  of  the  sum  II  is  valid  alto  for  II  . 

a  c 

The  total  inner  potential  is  obtained  by  adding  together  the  two  con¬ 
stituent:: 

n%  —  77f /7e.  (74) 

2.3  The  Boundary  Region.  ;  Between  0  And  -2c 

In  order  to  become  familiar  with  tfa-  sums  we  may  study  the  boundary 

region  in  which  both  expressions  II  and  should  be  valid.  Without  going 

U  *  J 

back  to  the  integrals  and  the  deformation  of  their  paths,  it  is  easily  shown 
on  the  series  themselves  that  both  expressions  are  identical  inside  the 
boundary  zone. 

To  this  end  we  assert 


nm-nt  =  n,,  (74a) 

and  have  to  prove  that  the  threefold  sum  on  the  right-hand  side  is  equal  to 
the  difference  of  the  double  s--ir-s. 
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It  is  easily  seen  that  the  third  summation,  over  n  on  the  right-hand 
side  has  its  origin  in  the  Fourier  development  of  the  parts  containing  u 
in  the  double  sums  of  the  left-hand  side.  In  compaxing  Eqs.  (71),  (72)  and 
(73)  one  sees  that  for  z  between  0  and  -2e  one  should  have  identically 
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(75) 


The  easiest  way  to  prove  this  relationship  is  to  prepare  the  sum  on  the 
righthand  side  in  the  following  way:  Start  with  the  development  of  e~ir:  in 
the  interval  0  . . .  -2c  which  is 
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Bring  the  factor  sin  vc  •  e  over  to  the  other  side;  there  remains  a 
development  with  coefficients 


The  coefficients  of  the  series  in  Eq.  (75)  are 
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Therefore  the  value  of  the  series  in  Eq.  (75)  is  the  difference  of  two 
exponential  functions: 
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The  righthand  side  of  this  equation,  completed  by  a  factor  i/c  is  after  some 
easy  manipulation  equal  to  the  lefthand  side  of  Eq.  (75).  This  shows  the 
identity  of  both  developments,  namely  Eq.  (71)  and  Eq.  (72)  +  Eq.  (73) 
within  the  boundary  region. 


3.  APPROXIMATIONS  FOR  OPTICAL  ft A\  ES 

Before  we  proceed  to  satisfy  the  dynamic  conditions  at  each  dipole  by 
means  of  the  incident  optical  wave,  we  will  attempt  to  obtain  a  survey  of 
the  nature  of  the  fields  represented  by  the  potentials.  There  are  two  cases 
to  distinguish: 

(1)  In  the  case  of  optical  waves,  the  quantify  a/A,  that  is,  the  half 
distance  of  the  lattice  points  measured  in  wavelengths,  is  a  small  number 
of  the  order  of  magnitude  1/100.  In  our  equations  this  means  that 

k  =  n/q  =  2VA  and  together  with  it  kQ,  a,(i, y  are  all  small  compared  to  1/a. 

(2)  If  an  X-ray  goes  through  the  lattice,  the  rigorous  computation 

based  on  the  assumption  Eq.  (61)  is  valid  as  in  case  Eq.  (58).  This  time, 

however,  ak  is  of  the  order  of  magnitude  of  10,  since  the  wavelength  of 

o 

the  usual  X-rays  (maximum  1A)  is  more  than  1000  times  smaller  than  the 
wavelength  of  light.  This  case  will  be  treated  in  Section  7. 

3.1  Optical  Waves 

3. 1. 1  THE  INNER  WAVE  II  . 

Q 

The  potential  II  is  already  known  from  Part  I.  It  is  sufficient  here  to 

fl  iks 

recall  that  it  contains  two  factors,  one  of  which  has  the  form  e  of  a  plane 

wave  which  proceeds  in  the  same  manner  as  the  dipole  phase  of  the  lattice. 

The  other  factor  is  a  threefold  sum  which  describes  in  the  manner  of 

Fourier  a  function  which  repeats  from  one  elementary  parallelepiped  to 

another.  The  optical  field  must  obviously  be  defined  by  an  averaging  process 

in  such  a  way  that  the  variations  of  this  rapidly  changing  constituent  are 
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eliminated.  For  this  purpose,  one  integrates  the  second  factor  over  a 
parallelepiped  thus  obtaining  an  average  potential 


nt  =  n™  = 


n  -Hal  +  fhj  +•  »it  1 
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(77) 


To  see  how  this  is  related  to  the  macroscopic  field  strength,  see  Section  5. 

3. 1. 2  THE  INNEB  WAVE  II  , 

c 

The  nature  of  the  terms  of  this  double  sum  is  determined  by  the  value 
of 


= + fkT- 


In  the  optical  case  only  the  term  (0, 0) 


(78) 
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can  be  real:  Ail  other  terras  are  definitely  imaginary.  All  terms  of  the 
sum  except  for  i  «  m  =  0  represent  inhomo-;  .meoue  waves,  whose  ampli¬ 
tudes  decrease  in  the  z  direction  according  to  the  factor 


elTim- 


The  higher  the  indices,  so  much  fester  *s  the  decrease,  if  we  have 

v0l  =  W— «*-(»/ *-fOs  ~ 
the  factor  becomes 


(79a) 
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If  we  take  the  case  of.  cubic  lattice,  where  a  =  b  =  c  it  thus  occurs  that  the 

amplitude  of  the  wave  through  the  first  dipole  layer  z  =  0  has  only 
_2r 

e  -  0.0018  of  its  amplitude  at  the  start  z  =  -2c  of  th-*  bcunJary  region. 
The  same  ratio  of  the  amplitudes  of  the  inhomogeneous  wave  holds  between 
the  first  and  second  dipole  layers. 

The  rapid  decrease  of  the  inhomogeneous  waves  towards  the  inside 
allows  one  to  assume  that  they  are  of  no  impoitance  for  the  essential 
process  of  reflection  and  refraction.  They  do  not  noticeably  influence  the 
dynamic  condition  of  the  lattice  even  in  the  outermost  layer  of  atoms. 
Besides  the  attenuation  towards  the  interior,  also  the  amplitude  itself  of 
the  inhomogeneous  waves  is  of  an  entirely  different  order  of  magnitude  than 
that  of  the  homogeneous  waves,  especially  the  homogeneous  wave  00  of  the 
potential  IIc.  For  the  term  00  the  factoi  l/i^  in  Eq.  (73)  ia,  because  of 
Eq.  (79)  of  the  order  of  magnitude  lA  .  For  the  inhomogeneous  waves  of 
the  higher  terms,  however,  it  is  only  of  magnitude  1  Eq.  (79?.) . 

From  all  this  it  is  justifiable  to  neglect  the  inhomogeneous  waves  for 
the  whole  inside,  including  the  atoms  in  the  boundary  and  to  take  into  con¬ 
sideration  for  the  dynamic  condition,  which  must  be  satisfied  at  each  atom 
only  the  potential  and  the  term  00  of  the  potential  II,. 

We  must  briefly  inspect  this  last  term.  We  have 


*oo  =  yv-«s-/s*  =  VV-**  +  r* 

*  — r*  sin !  3  > 


—  i0  cos  30 


where  v  =  k/k^  is  the  index  of  refraction,  and  (8)  and  are  respect¬ 
ively  the  angles  formed  with  the  z  axis  by  the  refracted  ray  (constructed 
according  to  the  law  of  refraction)  and  by  the  incident  ray: 


«n  3, 

sin  3 
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The  only  remaining  term  from  IIc  is,  when  written  out  in  full: 


Boo 


in  iei  +  i^y  +  ito  eo«  S»t  1  1 

2 ah*  1  _e-2Mkoe°«^-r)  ^  cos  8, 


(80) 


The  wave  normal  thus  lies  in  the  plane  formed  by  the  refracted  wave  and 
the  z  axis  (plane  of  incidence) .  The  phase  velocity  of  the  wave  will  be 
recognized  from 


\f«rTWTl -  k0  [See  Eq.  (79b; ] 

as  the  light  velocity  in  free  space.  Finally  the  angle  to  the  z  axis  is  the 
same  angle  (9^),  under  which  an  incident  wave  would  produce  the  internal 
wave  by  refraction. 

In  the  approximation  which  neglects  term£  of  higher  order  in  ck^  the 
factor 

_ l _ 

may  be  replaced  by 
_ i _ 1 

He  i*  cos  S,  -  r  (80a) 

_  _1 _ 1 _  _1 _ 1_  «n  3 

2te  k,  coe  S,  —  k  coe  3  2t'e  i,  «n  l3  —  8J 

3.1.2  THE  EXTERNAL  WAVES  Jl. 

a 

The  potential  Eq.  (71)  in  outer  space  is  very  similar  to  the  potential  IIC 
just  discussed.  Also  here  the  term  00  exceeds  the  amplitude  and  significance 
of  the  other  terms  and  one  can  restrict  oneself  to  that  one  term. 


Figure  8 
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4.  FULFILLMENT  OF  THE  DYNAMIC  CONDITION?, 

THE  FRESNEL  FORMULAS 

The  theory  of  the  unbounded  crystal  has  as  its  objective  to  show  that 
the  assumption  of  planes  of  equal  dipole  phase  leads  to  a  dynamically  pos¬ 
sible  state.  As  a  condition  of  this  self-consistency  it  was  found  that  the 
propagation  velocity  of  dipole  phase  and  with  it  that  of  its  accompanying 
light  wave  must  have  a  value  determined  from  the  dispersion  formula. 

For  the  bounded  crystal  the  field  by  which  a  dipole  is  excited  results 

from  t/v-o  potentials,  the  potentials  II  and  II  ,  from  which  however  there 

q  c 

is  to  be  subtracted  the  part  which  is  emitted  by  the  excited  dipole  itself. 


Since  II 
Q 


everywhere,  even  for  the  outermost  boundary  atoms,  is  identical 


with  the  potential  in  the  infinite  crystal,  the  same  field  of  excitation,  and 


therefore  the  same  mode  of  dipole  oscillation  is  generated  in  the  bounded 
as  in  the  unbounded  crystal  only  if  Ec  can  be  suppressed.  As  the  dynamic 
conditi  -u  for  the  feasibility  of  the  assumption  of  planes  of  equal  dipole 
phase  we  retain  here  as  there  only:  that  one  must  choose  the  correct  phase 
velocity  q.  In  the  optical  case  the  potential  II  reduces  to  a  single  homo- 
geneous  wave  that  geometrically  coincides  with  the  incident  wave  extended 


into  the  inside  of  the  crystal.  Therefore  the  amplitude  and  phase  of  the 
dipole  wave  can  easily  be  so  adjusted,  that  C  and  the  incident  wave  cancel. 

In  order  to  accomplish  the  transition  from  the  potential  to  the  field 
strengths,  we  take  some  considerations  beforehand  which  will  be  system¬ 
atically  discussed  in  Section  5.  In  doing  so  we  limit  ourselves  here  to  the 
case  of  a  cubic  crystal  (a  =  b  =  c,  wave  velocity  in  all  directions  of  equal 
value,  dipole  oscillation  at  right  angles  to  the  propagation  direction  s; 
compare  dispersion  theory).  All  of  the  potentials  when  completely  written 
out  lave  the  form 


eaein* , 

where  x  is  either  k  or  1^  and  ea  =  p  signifies  the  moment  of  the  dipole. 
According  to  Part  I,  Section  2,  Eq.  (4),  the  field  strength  is  obtained  from 


the  potential  by  tbe  operation  curl  curl.  We  assert  here  and  prove  it  in 
Section  5,  that  this  operation  in  the  general  case  yields  tbe  dielectric  dis¬ 
placement  3> .  Only  in  the  case  of  empty  space,  £  =  1,  does  it  likewise 
yield  the  electric  field  strength  £,  which  then  is  equal  to  X . 

Take  as  tbe  propagation  direction  tbe  x  axis;  we  obtain  from  the  above 
potential 


=  0, 

Or  z  or  1 

c  0  -  _  .  0* 


_  £  q  1  _ A 

by  x  dr  *  Si1 


"dy  *  a7 

T7  a*7r*m*  -  *•£* 


0t~rei«  =^0^-. 


This  means  generally  expressed:  In  converting  to  tbe  field,  only  the 
component  of  the  amplitude  at  right  angles  to  the  wave  normal  produces  a 
centribi£Ion.  Tbe  field  of  tbe  vector  3)  differs  from  the  field  of  the 
transverse  pari  cf  tbe  vector  potential  by  the  additional  amplitude  factor 
x*  (  =  k2  or  kg2).  We  have  thus  to  distinguish  two  cases:  one  in  which  the 
dipole  oscillations  occur  at  right  angles  to  the  plane  of  incidence,  and  the 
second  in  which  the  oscillations  occur  in  tbe  plane  of  incidence.  fo  tbe  first 
case  tbe  full  amplitude  a  is  valid  for  all  three  waves;  in  the  second  case 
tbe  amplitude  is  reduced  by  a  cosine  factor. 

From  the  potentials  Eqs.  (77),  (80),  and  (81)  we  obtain  according  to 
this,  the  following  field  strengths  for  tbe  refracted  wave,  tbe  boundary 
wave,  and  the  reflected  wave  when  tbe  oscillations  occur  in  tbe  plane  of 
incidence: 

L  Refracted  wave: 


35  =  — a— 

»  * 


t  k,*-k* 


*  The  abbreviation  tg  has  not  been  changed  in  these  equations  —  it  is 
obviously  tan. 

**  The  following  formulas  are  useful  for  the  conversion: 

1  -  **  =  co*3  3,  -  r1  co*3  3 

*inf3  -  3*}  *j2{3  —  3g) 
sia*  3 

eo*  3,  -  *  cc*  3  _  ea (3-jJ 
co*  3,  -  r  eo*  3  cn  {3  3*) 


Further: 
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snH  the  amp!  finite  of  the  reflected  wave  as 

?3>J  _  gafS+y  gaCS-St)  tg(3»-3) 

*  c©»{5-&*}  *is  {3+3*3  “  Jg(S  +  3^  *  *,**  (83b) 


Finally,  for  tie  case  in  which  the  oscillations  are  polarized  at  right 
angles  to  the  plane  of  incidence,  the  refracted  rave  would  be  mult'-plied  by 
the  factor  cos  (8  —  8*)  and  the  reflected  wave  would  be  multiplied  by 


co«  (3  —  3,) 

«•(*  +  *») 

The  formulas  for  oscillation  at  right  angles  to  the  plane  of  incidence  are 
thus: 


*1 


2- 


1^7 


>(3- V 


rfaS 


►coa  8,  *  2**  co«  S#- 


»{»  +  V 


(84a) 


J5.  «a(S  +  S.T 


(84b) 


These  formulas  are  in  agreement  with  the  Fresnel  formulas  for  the  vector 
3) ,  the  dielectric  displacement  (compare  the  formula  for  %  in  Drude  Lehrb. 
d.  Optik,  1.  Edition,  p.  258-262). 

The  preceding  treatment  of  the  problems  of  reflection  and  refraction 
has  been  developed  strictly  within  foe  framework  of  the  theory  of  electrons. 
&  confirms  foe  expectations  of  the  theory  of  dispersion  regarding  the  in¬ 
fluence  of  the  boundary.  The  feet  that  an  internal  boundary  wave  appears  in 
the  electron  theory  is  the  point  which  allows  a  deeper  insight  into  the  prob¬ 
lem  than  the  derivation  according  to  Maxwell  on  the  basis  of  boundary  con¬ 
ditions.  Even  though  these  boundary  conditions  were  established  from  the 


*,  **See  footnote  on  page  68. 
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concepts  of  electron  theory,  the  very  characteristic  property  of  tbt 
boundary  lies  concealed  In  the  boundary  conditions  and  it  does  not  appear 
as  distinctly  as  above:  why  the  electrons  are  bound  to  oscillate  in  the 
described  manner  in  order  that  they  harmonize  with  the  total  field  of  excita¬ 
tion. 

The  difference  of  the  Fresnel  formulas  in  the  two  cases  of  polarization 
is  clarified  here  in  a  very  direct  way  as  being  due  to  the  component  forma¬ 
tion  of  the  dipole  amplitudes.  A  particular  case  of  this  is  the  familiar  ex¬ 
planation  of  the  Brewster  angle,  as  that  angle  of  incidence  for  which  the 
oscillations  o?  electrons  in  the  body  (so  far  as  they  occur  in  the  plane  of 
incidence)  have  the  direction  of  the  reflected  ray  and  are  thus  optically 
ineffective  in  producing  its  amplitude. 


5.  OBSERVATIONS  ON  THF  TRANSITION  FROM  THE  MICROSCOPIC  FIELD 
TO  THE  VECTORS  g  AND  %  OF  THE  MAX7ELL  THEORY 

The  Fresnel  formulas  were  derived  in  the  preceding  paragraphs  by 
replacing  the  potential  11^  by  its  average  value  over  an  elementary 
parallelepiped;  by  farther  performing  the  operation  curl  curl  on  this  average 
value,  and  by  claiming  that  the  macroscopic  field  of  the  vector  3)  may  be 
obtained  by  this  method  from  the  microscopic  field. 

The  reason  for  this  procedure  is  analyzed  in  the  following. 

The  vector  potential  was  chosen  in  view  of  an  easy  description  of  the 
field  of  a  dipole,  fa  so  doing,  in  particular  by  setting  g  =  curl  curl  % 

(Part  I,  Section  2,  Eq.  (4)),  it  was  assumed  that  the  field  propagates  in 
empty  space  («  =  1).  The  potential  e^^/R  satisfies  the  differential 
equation  for  vacuum 

=  (85a) 

By  this  rule  also  in  Part  I,  Section  7,  the  excitation  field  was  obtained 
from  the  exciting  potential. 


. . .  '  *'•  ‘ . . . *w<iiiMMii)ii(iuiuuMUu)DUUIUUbt(l|QtlUliltUIIII(UIMHnHltHliCII. 
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By  averaging  the  total  potential  a/7  to  smooth  the  fast  fluctuations, 
we  find  that  the  optical  potential  now  satisfies  the  equation 

J$-rA*$=0  (55b) 

2 

k  is  the  constant  occurring  for  the  propagation  in  a  medium  with  the 
dielectric  constant  s  (=  u2).  But  whenever  a  vector  potential  satisfies 
this  equation,  then  the  operation  curl  curl  produces  fro'n  it  the  vector 
the  dielectric  displacement. 

This  follows  from  the  definition  of  the  vector  potential.  It  is  intro¬ 
duced  for  a  body  with  the  dielectric  constant  «  so  as  to  lead  to  a  solution 
of  Maxwell’s  equations 


Ve  —  —  curl 

S; 

(S6a) 

0. 

II 

(86b) 

div  $  =  0 

(86c) 

div?)  —  0 

(86d) 

(86e) 

One  determines  in  view  of  Eq.  (86c)  curl  8  such  that 
curl  8  =  $ 

the  second  equation  then  makes 
=  curl  curi  8, 

whereby  also  Eq.  (86d)  is  satisfied. 

One  can  now  dispose  of  div  8  by  putting 


div  8  =  0, 


thus  Eq.  (86a)  in  connection  with  Eq.  (86e)  requires 


or  JS  +  A**-0.  (85b') 

Let  it  be  noted  immediately  that  for  an  orthorhombic  crystal  the  relation¬ 
ship  between  ©  and  3)  is  not  given  by  Eq.  (86e)  but  by  the  three  equations 

3>,  =*  e,®,;  3).  =*  «,€,  (86e!) 

and  that  thus  the  three  components  of  31  satisfy  different  equations 

Eq.(85b),  in  which,  instead  of  the  uniform  e  are  substituted  e  ,  ev,  t  for 

X  j  z 

the  three  components  But  here  as  in  the  isotropic  body 

(cubic  crystal)  it  is  always  true  that 

curl  curl  =  5).  (85c) 

So  much  for  the  reasons  that  the  method  used  provides  the  Fresnel 
formulas  for  the  vector  35  and  not  for  ®. .  However,  a  remark  may  be 
added  on  the  transition  from  the  differ  atial  equations  Eqs.  (85a)  to  (85b), 
which  is  connected  with  the  definition  of  the  values  ©  and  3)  and  which 
throw's  an  interesting  light  on  the  logical  aspect  of  the  fir  st  part  of  the 
theory  of  dispersion. 

The  total  potential  satisfies  Eq.  (85a)  before  averaging:  Ea.  (85b)  after 
averaging.  The  integration  over  the  elementary  parallelepiped  can  pro¬ 
duce  a  change  in  the  differential  equation  only  because  it  includes  an  in¬ 
tegration  over  the  pole.  For  the  integration  of  ,  which  would  become 
infinite  as  l/R,  the  pole  is  insignificant.  However,  if  is  integrated 
over  the  pole,  the  near  vicinity  of  the  pole  contributes,  as  is  well  known, 
a  term  4 7tp,  where  p  ~  e*  a  signifies  the  moment  of  the  dipole,  while 
the  remainder  of  the  integration  space  produces  A  where  $  denotes 
the  average  potential.  The  average  of  the  Eq.  (85a)  over  the  volume  v  thus 
gives: 
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+  (85d) 

V 

If  now  the  average  of  the  potential  is  the  potential  of  the  Maxwellian 
quantity  S5  then  it  must  also  satisfy  the  propagation  equation  for  %  . 

That  is,  we  must  also  have 

A%  +  «  0.  (85e) 

From  these  two  equations,  the  average  value  of  the  potential  itself  follows: 

(V -  A*) f5  =  ~^rP  -  — 

or 

(85f) 

But  this  is  exactly  the  result  from  the  theory  of  disoersion  [  Part  I, 
Section  6,  Eq.  (31)]  which  was  again  derived  in  Part  n  by  averaging  ovei 
the  sum  of  the  dipole  fields  [Eq.  (77)]. 

One  sees:  If  it  Is  assumed  that  the  average  value  of  the  potential  is 
the  potential  of  the  optical  field  and  therefore  satisfies  Eq.  (85e),  then  it  is 
cot  necessary  first  to  calculate  the  exact  fields  produced  by  the  dipole 
oscillations.  Eather,  one  can  infer  the  Maxwell  field  quantities  directly 
from  the  assumption  made  concerning  the  dipole  oscillations.  Whether  the 
assumed  state  of  dipole  oscillations  is  dynamically  possible,  one  is,  how¬ 
ever,  not  in  a  position  to  decide  without  going  into  the  details  of  the  field. 
But  only  a  knowledge  of  the  fields  in  the  immediate  vicinity  of  the  dipoles 
(the  excitation  fields)  wouid  be  required  for  finding  the  answer. 

The  logical  place  of  our  computation  of  the  total  field  from  the  as¬ 
sumed  dipole  oscillations  is,  in  Part  I,  the  theory  of  dispersion,  an 
existence  proof.  Without  using  the  assumittions  necessary  for  the  deriva¬ 
tion  of  Eq.  (85f),  it  is  shown  that  there  exists  an  average  value  of  the 
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potential  which  satisfies  Eq.  (85b)  anti  is  therefore  the  potential  of  the 
optical  field.  In  the  theory  of  reflection  and  refraction  on  the  other  hand 
the  computation  of  the  total  field  serves  the  purpose  of  resolving  the 
dynamics  of  the  problem.  It  could  therefore  not  be  omitted,  even  if 
relationship  Eq.  (85f)  were  used. 

In  conclusion  we  compare  our  definition  with  that  used  by  H.  A. 
Lorentz.  ^  The  latter  is,  written  in  our  units 

®  ®  +  4»  p, 


or 


$  =  •  4»p  *  -  4*  •  (85g) 

The  same  relationship  between  the  dielectric  displacement  and  polarization 
we  also  obtain  from  the  relationships  Eqs.  (85c)  and  (85f)  in  all  cases  in 
which  the  assumption  of  the  average  polarization  p  is  descriptive  of  a  state 
which  is  consistent  with  the  Maxwell  Eqs.  (86a)  to  (86d),  that  is,  whenever 

Jp+^pssO;  divp  =  0 

We  have  then  from  the  definition  Eq.  (85c) 


3)  **  curl  curl  % 


—  4  *  curl  curl  jj 


in  agreement  with  Lorentz. 


(1)  Enzyklop,  d.  math.  Wiss.  V,  14  (XXI),  (XXV),  (XXX3)  and 
No.  30. 
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The  definition  contained  in  Eq.  (85c)  thus  leads  to  the  familiar  vector  of 
dielectric  excitation  or  displacement  for  all  states  of  dipole  motion  which 
permit  comparison  with  a  Maxwellian  field. 


6.  WAVE  EXIT,  TOTAL  REFLECTION 


In  order  to  study  the  exit  of  a  plane  wave  from  the  body,  which  as 
before  fills  the  upper  half  space,  only  small  changes  in  the  final  formulas 
are  required.  This  case  differs  from  the  preceding  one  only  by  the  reversal 
of  the  wave  direction,  that  is,  the  values  a,  (3,y  have  to  be  replaced  by  their 
negative  values.  We  make  this  change,  and  the  principal  terms  of  the  three 
potentials  Eqs.  (77),  (80)  and  (81)  then  read: 


-iax-ity  -i*mi 

FT  ,n  1  1 

•  ~  2abe  !  _  ’ 


(87a) 


rt  _  tn  _ } _ _  _i_ 

e~2abe  i  _  e-2»c(**»  +  r)  rm  1 


(87b) 


-iax-./y-iy*  , 

IT  71  1 

11  n  "  2abc  *,*-** 


(87c) 


Here  has  the  old  value: 


*00  «  VV-«*  *o  Vl-^’S .  (87d) 

In  nq  we  recognize  again  the  old  wave  of  velocity  q,  now,  however, 
running  towards  the  boundary.  Its  amplitude  is  unchanged. 

JIa  as  well  as  IIc  have  the  phase  velocity  c  and  have  their  angle  with 
reference  to  the  z  axis  preserved,  but  the  angles  to  the  x  and  y  axes 
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changed,  so  that  the  negative  axis  directions  replace  the  positive  ones. 
Also  the  previously  given  amplitudes  are  changed.  We  substitute  as  above 
[Section  3  Eq.  (80a)] : 

- I -  by - i - 1 — , 

1  _  e  2«'c(*«.-r  2  tcrm-j 

(87’) 

_ J _  by  J-  _L-  . 

I  _e  *2»c(rw  +  r»  J  2*c>-„+r 

(87”) 

The  pictures  of  the  waves  which  result  from  incidence  and  exit  are 
compared  with  each  other  in  Figure  9. 


Figure  9 

One  now  sees  easily  how  to  fulfill  the  dynamic  requirements.  Tie 
internal  wave  of  velocity  c  has  to  be  removed,  since  ;t  would  disturb  tie 
electrons  whose  oscillations  are  balanced  only  against  the  wave  of  velocity 
q.  This  is  accomplished  by  considering  the  dipole  oscillations  coupled  to 
tie  internally  reflected  wave  q’  which  would  produce  an  internal  boundary 
wave  c’  of  the  same  direction  as  c.  (Figure  10).  Its  amplitude  is  deter¬ 
mined  by  the  requirement  that  c  and  c’  cancel  each  other.  This  gives  at 
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and  hence  from  Eqs.  (87a)  and  (87'),  and  Eqs.  (77),  (81),  and  (81a) 


n  =  a  4-  a  ... 


1  _J - L(*  *  _  k- )  -  —  (*0l- A*). 

*  - 1  'oo  K  0  ("«.  +  7> 


In  going  over  to  the  fields,  q  and  q’  have  to  be  multiplied  with  the  factor  k  , 

2 

the  wave  a"  in  outside  space  however  by  kg  .  If  we  continue  to  set  the  am¬ 
plitude  of  the  field  strength  of  the  wave  q  equal  to  1,  then  the  amplitude  a” 
contains  the  factor  k^  /k  =  l/v  .  3y  utilizing  the  formulas  mentioned  in 
the  footnote  section  4,  page  68,  the  amplitudes  of  the  internally  reflected 
and  the  transmitted  wave,  for  the  case  that  the  oscillation  is  at  right  angles 
to  the  plane  of  incidence,  become 


»in  (3-3J  _2_  cos  (3  rin 

sin(3+3,j'  a  »*  sin (S  +3^ 


(89) 


These  are  the  same  values  as  are  obtained  from  the  Fresnel  formulas 
Eqs.  (84a)  and  (84b)  when  the  angles  (8)  and  (8,)  are  interchanged  and 
v  is  reolaced  by  l/v  corresponding  to  the  progression  from  an  optically 
denser  to  an  optically  thinner  medium. 

Fulfillment  of  the  dynamic  conditions  has  thus  provided  the  correct 
amplitudes  also  in  the  case  of  wave  exit.  It  was  tacitly  assumed  for  the 
last  conversions  that 


Agf'l  —  r*  sin*  3 


may  be  replaced  by  cos  85.  This  implies  that  we  are  not  dealing 
with  the  case  of  total  reflection. 

Let  us  assume  now,  that  the  angle  of  incidence  (8)  of  the  internal 
wave  with  the  boundary  surface  is  so  large  that  total  reflection  occurs. 
i/qq  is  then  imaginary  and  the  two  boundary  waves  associates  with  the  q 
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wave  in  inner  and  outer  space  are  no  longer  homogeneous,  but  inhomo¬ 
geneous.  They  agree  in  this  regard  with  the  higher  terms  of  the  series 
Eqs.  (71)  and  (73)  for  the  potentials  II  and  II  whose  principal  (00)  terms 
they  were.  But  compared  with  these  higher  terms,  for  instance  the  term 
(01),  they  show  a  quantitative  difference  in  the  degree  of  inhomogeneity. 

While 

*09  =  W  —  a* -{z  I  ft- [if 

has  an  imaginary  value  of  the  order  of  magnitude  of  l/b , 

—  ik6  \  r-  sin1  3  —  1  ~  i  ^ 

is  only  of  the  order  of  magnitude  k^,  flat  is,  1/K  The  inhomogeneous 
wave  (00)  thus  reaches  with  noticeable  strength  as  far  down  in  depth  as  the 
value  of  the  wavelength,  while  the  higher  order  waves  do  not  reach  beyond 
the  first  atomic  layer.  Also,  as  before,  since  the  factor  l/y^  appears  in 
the  absolute  vain  ,  of  the  amplitude,  the  higher  inhomogeneous  waves  are 
much  weaker  than  the  principal  wave  (00).  Physically  speaking,  the  higher 
inhomogeneous  terms  describe  the  perturbation  at  the  surface,  which  is 
caused  by  the  discontinuous  nature  of  the  surface.  The  (0C)  wave,  on  the 
other  hand,  is  produced  by  the  events  in  the  interior  and  is  independent  of 
the  smallness  or  coarseness  of  the  surface  lattice  of  atoms. 

The  dynamic  conditions  remain  the  same  for  the  total  reflection  as  for 
the  previously  discussed  exit  of  a  wave  from  the  medium.  The  internal  in¬ 
homogeneous  wave  c  must  be  removed  by  considering  the  internally  re¬ 
flected  wave.  Inside  the  medium,  the  two  inhomogeneous  waves  c  and  c' 
are  superimposed.  For  the  reflected  wave  in  the  interior  we  obtain  as 
above  Eq.  (38b)  the  amplitude 

rm  —  r  k  cos  3  —it,  W1  sin*  3  —  1 

—  * —  ITS.  ■  ■ 

**•  +  7  k  COS  3  -i-  it.  V**  sin*  3—  3 
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Its  absolute  value,  which  is  the  important  quantity,  is  ?  since  numerator 
and  denominator  are  complex  conjugates. 

The  potential  of  the  inhomogeneous  wave  in  outer  space  has  at  the 
boundary  the  amplitude  [  compare  Eq.  (88c)  ] 

X  J_fA  _  *m-r)  _ _ _ 

'-~7  Irm  +  f))  (*Vr*KB*S-l  +  t  co*  S) 

The  absolute  value  of  this  is 


2  cos  s 


2  2  2 

For  the  field  of  the  electric  vector,  we  have  to  add  the  factor  /k  =  \/v 
for  the  same  reason  as  previously.  That  the  result  coincides  with  the  one 
obtained  from  Maxwell's  theory,  follows  from  the  fact  that  it  can  be  obtained 
here  as  above  from  the  /resnel  formulas  by  replacing  cos  8,  by 
t  V  Vsm*  8  -  i  . 

TRANSITION  TO  THE  CASE  OF  X-RAYS  * 

The  series  Eqs.  (71),  (72)  and  (73)  for  the  field  inside  and  outside  the 
crystal  are  valid  in  all  strictness  on  the  basis  of  the  assumptions.  The  dis¬ 
cussion  of  the  pervious  paragraphs  indicaies  that  for  the  case  of  optical 
waves  (ak  « 1)  the  three  series  reduce  essentially  to  their  principal  terms, 
in  which  all  indices  are  zero. 

it  is  otherwise  when  x-rays  penetrate  the  lattice.  **  We  consider  the 
field  in  outer  space,  whose  potential  is  the  series  Eq.  (71)  (The  same 
consideration  is  associated  also  with  the  exit  of  the  ^-rays.)  There  exists 

*  This  case  I  have  already  treated  on  the  basis  of  formula  (7  2)  in 
Phys.  Zeitschr.  14:465,  1913. 

-  VA 

**Tbe  lattice  is  saturally  thougiS  of  as  being  perfectly  regular,  hi 
particular  the  attenuating  thermal  motion  is  completely  neglected. 
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noy  a  whole  series  of  terms  for  which 

is  real,  and  which  therefore  describe  homogeneous  plane  waves.  Among 
these  waves  those  are  dominant  for  which  the  denominator  of  Eq,  (71) 
vanishes.  The  condition  for  its  first  factor  to  vanish  is  that 

c  \Tlm-1ry)=n7t 

where  n  is  a  positive  or  negative  integer.  Restated,  this  condition  becomes 
\l  zja  —  a?  +  (»  a/A  —  ft?  +  {*  zjc  —  7?  =  kft-  (30) 

Note  that  here  (certainly  to  a  high  degree  of  approximation) 


a*  -f  /?*  -  y*  =  V  (=  **)- 

Equation  (90)  leads  directly  to  fee  sphere  construction  for  the  deter¬ 
mination  of  the  direction  of  the  interference  maxima:  It  is,  in  fact,  identical 

(2) 

with  the  Lane  conditions'  for  the  occurrence  of  an  interference  maximum. 

That  the  amplitude  of  the  waves  appears  to  be  infinitely  large  is  due  to 
the  cooperation  of  the  infinite  number  of  dipoles  and  tie  neglect  of  absorption. 

Observe  moreover,  that  in  tie  interior  of  the  crystal  according  to  the 
,  potential  0Q  likewise  each  of  the  wave.-:  {/ ,  m,  n)  attains  a  maximum  strength 
for  which  Eq.  (90)  is  satisfied.  The  field  in  the  interior  consists  essen¬ 
tially  of  the  interference-favored  -waves,  *  s  cos  of  these,  and  not  feather 


(1)  F.  P.  Ewald,  loc.  cat.  or  EnrvkJon.  d.  math.  w:ss. ,  Article  Lane, 
V.  24,  No.  48. 

(2)  H.  v.  Lane,  Ann,  d.  Pays.  41;  STL  2nd  Enzrklon.  d.  math.  Wiss. 
Article  Laae  V.  24,  No.  46. 
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distinguished  from  the  others,  there  exists  the  (000)  wave,  which  in  the 
optical  case  was  the  only  dominant  one,  namely  the  refracted  wave. 

One  apparently  obtains  in  Eq.  (71),  an  additional  mode  of  infinite 
amplitude,  namely  when  the  second  denominator  vanishes: 

r,m  =  0. 


ft  follows,  however,  from  the  derivation  of  formula  Eq.  (71)  that 
application  of  Eq.  (71)  to  this  case  is  not  permitted,  (the  pole  is  in  this 
case  of  second  order). 

It  is  instructive  to  consider  the  transition  from  the  case  of  x-rays  to 
the  case  of  optics.  To  that  effect,  we  let  ak  gradually  decrease  —  we  can 
think  of  the  crystal  lattice  shrinking  while  we  held  the  wavelength  constant. 
There  are  thee  a  diminishing  number  of  index  pairs  i,  m,  which  leave 
real.  Keeping  i,  m  constant,  vIm  becomes  ever  smaller.  This  means  that 
the  fan  of  interference  rays  which  is  present  for  the  shorter  x-rays, 
spreads  out.  A  moment  comes  where  the  rays  which  are  diffracted  the 
most  hare  =  0.  From  there  on  for  these  rays  would  become  and 
remain  imaginary.  As  soon  as  the  wave  directions  approach  the  crystal 
surface  the  waves  become  inhomogeneous.  The  spreading  of  the  pencil  of 
interference  rays  can,  of  course,  be  understood  by  the  increasing  dispersive 
power  of  the  shrinking  lattice  which  spreads  the  spectra  apart.  Only  the  CO 
ray  (indicated  thick  in  Figure  11)  does  not  charge  its  position  in  this 


Figure  11 
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procet  s.  It  is  this  *“iy  that  finally  remains  as  the  only  one  of  the  reflected 
spectra  that  survives  in  the  optical  domain. 
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Postscript  1970 
by  P.  P.  Ewald 


DEVELOPMENTS  ARISING  GUT  OF  PARTS  I  AND  \i 
ON  THE  FOUNDATIONS  OF  CRYSTAL  OPTICS 

Before  the  end  of  the  summer  semester  1910  at  the  University  of 
Munich,  I  asked  my  revered  teacher.  Professor  Arnold  Soxnmerfeld,  for 
permission  to  start  on  thesis  work  under  his  guidance,  He  suggested  for 
subjects  some  ten  different  applications  of  Maxwell's  equations  to  boundary 
value  problem 8,  such  as  radio  wave  propagation,  or  the  self-inductance  of 
coils  of  various  shapes.  At  the  bottom  of  the  list  stood  the  problem:  Does 
anisotropic  arrangement  of  resonators  produce  double  refracticn  ?  This 
was  the  only  problem  that  ended  with  a  question  mark,  and  the  Professor 
explained  that  he  himself  saw  no  clear  way  of  tackling  the  problem  and  that 
it  would  be  much  more  risky,  from  my  point  of  view,  of  choosing  this 
problem  rather  than  one  of  those  where  he  could  foretell  what  the  solution 
would  require  in  the  way  of  mathematics.  He  explained  that  crystallo- 
graphers  had  long  since  developed  elaborate  theories  of  crystal  structure 
but  that  these  remained  purely  geometrical  speculations  as  long  as  no 
quantitative  consequences  had  been  drawn  *•  and  that  if  I  could  find  an  answer, 
it  might  help  to  corroborate  or  to  refute  the  ideas  on  crystal  structure. 
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Of  all  the  subjects  he  proposed,  only  this  one  evoked  any  enthusiasm 
in  me,  and  so,  after  thinking  it  over,  I  returned  to  him  with  the  decision 
that  I  would  like  to  work  on  it.  Professor  Sommerfeld  gave  me  the  re¬ 
print  of  a  paper  by  M.  Planck  on  the  theory  cf  dispersion  -  and  with  that  I 
started  on  a  solitary  hike  up  the  Rhine  valley  in  Switzerland  for  the  next 
weeks  of  summer  vacation. 

I  was  not  much  wiser  when  I  returned  to  Sommerfeld  in  the  fall.  Not 
only  had  I  not  progressed  with  my  problem,  but  I  even  found  great  difficulty 
in  understanding  what  Planck  was  doing,  and  Planck  had  always  been  praised 
to  me  as  a  master  of  clarity.  There  seemed  to  me  to  be  some  confusion  in 
his  calculation  of  the  field  created  by  the  resonators  (dipoles)  in  the  medium 
and  the  "incident  field"  which  seemed  to  exist  in  addition. 

Professor  Sommerfeld,  meanwhile  ,  lad  also  given  some  thought  to  the 
problem.  Using  the  Hertz  potential,  he  lad  tried  to  sum  up  the  fields  of 
the  spherical  waves  issuing  from  the  dipoles  at  the  nodal  points  of  the 
orthorhombic  iattiee  and  had  recognised  that  the  main  obstacle  in  this 
summation  was  the  fact  that  the  dipole  coordinates  X,  Y.  Z.  entered  in  the 
Hertz  potential  only  in  the  irrational  expression 

R  =  ^(X-x)2  +  (Y-y)2  4-  (Z-z)2 

where  (x,  y,  z)  denote  the  field  point.  He  then  shewed  me  a  representation 
of  by  an  integral  containing  Bessel  functions  of  r,  the  two-dimensional 

analogue  to  R,  and  an  exponential  factor  in  which  the  exponent  was  linear 
in  the  third  coordinate  (Z  -  z>.  The  summation  over  the  third  coordinate 
could  then  be  readily  performed  under  the  integral  sign,  for  it  was  nothing 
but  the  sum  of  a  geometrical  series.  Cculd  we  net,  Sommerfeld  suggested, 
bring  the  other  two  coordinates  into  a  similar  linear  form  ?  For  the  next 
few  months,  I  triad,  in  vain  to  achieve  this  -  then,  suddenly.  I  saw  the 
light  and  obtained  the  desired  representation  in  the  form  of  the  triple  in¬ 
tegral  (Part  I,  Section  5,  Eq.  (23)].  When  I  arranged  to  show  Sommerfeld 
my  recent  advance,  he  greeted  me  with  the  words:  '1  too  can  show  you 
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something  that  should  be  helpful1'  -  anc  he  developed  the  very  same  integral 
which  he  had  derived  in  a  less  cumbersome  way. 

The  path  now  iay  open  for  finding  the  total  potential  by  means  of  the 
residue  method  and  I  rapidly  arrived  at  the  expression  Part  I,  Section  6, 

Eq.  (30), 

My  happiness  over  this  first  tangible  result  was  not  to  last  long:  I 
found  it  impossible  to  subtract  from  this  total  potential  that  part  that  was 
contributed  by  the  test  dipole.  Since  the  oscillations  of  the  test  dipole  are 
caused  only  by  the  radiation  coming  from  the  other  dipoles,  this  subtrac¬ 
tion  appeared  essential  for  the  determination  of  the  "field  of  excitation". 

The  impossibility  of  a  direct  subtraction  comes  from  the  fact  that  in  the 
expression  Part  I,  Section  6,  Eq.  (30),  all  individual  dipole  fields  are 
jumbled  together  and  recast  into  a  sum  cf  norJ.ocamed  plane  waves  in 
which  individual  contributions  can  no  longer  be  recogrdzed.  (I  was  rot 
aware,  at  the  time,  that  what  I  had  obtained  was  a  Fourier  development  of 
the  iota,  field.) 

Success  in  the  removal  of  the  field  of  the  test  dipole  came  through 
Somme rfeld's  assistant,  P.  Debye.  At  a  skiing  holiday  in  Mittenwald, 
Easter  1911,  Sommerfeld  showed  Mm  our  quandary.  With  one  glance  at 
the  expression  of  the  total  potential,  Debye  said:  "This  is  quite  simple;  you 
have  to  use  Biemann’s  method  of  bringing  the  denominator  of  the  sum  term 
into  the  exponent  of  an  exponential  function  by  introducing  a  new  integration. 
The  integrand  then  becomes  a  theta-function  and  you  can  apply  the  trans¬ 
formation  theorem  of  the  theta  functions  to  it.”  The  whole  "consultation" 
probably  lasted  no  more  than  a  q’iarter  of  an  hour.  Debye  was  not  only  an 
avid  reader  of  classical  physics,  but  also  had  an  extraordinary  power  of 
seeing  through  mathematical  formalism.  He  did  not  know'  at  the  time  in 
which  of  Riemann’s  papers  this  method  occurred  -  nor  did  I  ever  find  out. 

When  I  Mid  finally  understood  Debye’s  advice,  I  transformed  the  total 
potential  in  the  way  he  had  suggested.  This  is  a  rather  round-about  way  of 
obtaining  the  theta  functions  and  in  the  publication  I  followed  a  more  direct 
way  by  introducing  the  simple-integral  representation  of  ^  [  Part  1, 

Section  5,  Eq.  (18)3,  wMch  leads  directly  to  theta  functions. 


A4 


For  an  orthorhombic  lattice,  the  theta  function  in  (X,  Y,  Z)  is  the 
product  of  three  such  functions  of  one  variable  each,  and  the  transformation 
is  easily  written  down.  The  integrands  in  the  two  forms  thus  obtained  are 
equivalent  at  any  point  of  the  integration  over  e .  In  fact,  one  ts  the  Fourier 
development  of  the  other.  It  is,  therefore,  permissible  to  use  one  integrand 
up  to  an  arbitrary  point  E  of  the  whole  integral  path,  and  then  change  to  the 
other  form  for  the  remainder  of  the  integration.  In  the  integral  from  E  to  » 
individual  contributions  to  the  potential  are  still  recognizable,  and  the  con¬ 
tribution  of  the  test  dipole  can  be  omitted;  whereas  in  the  integration  from 
0  to  E,  the  Fourier-transformed  integrand  is  used  in  which  individual  con¬ 
tributions  are  blotted  out.  By  suitably  choosing  the  transition  point  E,  one 
arrives  at  having  rapid  convergence  of  the  sums  in  both  integrals.  A  par¬ 
ticular  advantage  lies  in  the  fact  that  by  shifting  E,  rapid  convergence  can 
be  shifted  from  one  sum  to  the  other.  Since  the  sum  of  the  two  part  integrals 
is  independent  of  E,  an  easy  check  on  the  correctness  oi  tie  numerical  work 
is  afforded  by  using  two  different  values  E. 

I  have  later  extended  this  method  of  calculating  lattice  potentials  and 
Madelung  constants  from  orthorhombic  to  general  lattices.  There,  the 
theta  function  of  the  three  space  coordinates  no  longer  splits  into  the  pro¬ 
duct  3(x)  •  0(y)  •  6{ z>.  The  transformation  property  of  the  space-theta  was 
contained  in  the  main  text  on  theta  function,  * **  but  in  a  horrible  form.  By 
using  the  concept  of  the  reciprocal  lattice,  I  could  re-state  it  in  a  way  akin 
to  the  tools  of  the  mathematical  physicist.  Since  then,  the  theta  function 
method  has  been  the  favored  one  for  calculating  lattice  energies  and  po¬ 
tentials.  The  neat  way  in  which  the  formal  splitting  cf  the  integral  for  the 
potential  produces  two  rapidly  convergent  parts  was  interpreted  physically 
in  a  later  paper  of  mine.  *** 


*  Ann.  d.  Phys.  Lpz.  (1921}  64:253-287. 

**  A.  Krazer  Lehrbuch  der  Thetafunktionen,  Teubner,  Lpz.  1903. 

Nachr,  Ges.  d.  Wiss.  zu  Gottingen  (1938)  Math-phys.  Kl.  3:55-64. 
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For  bringing  my  thesis  to  a  conclusion,  two  points  had  to  be  discussed: 

(1)  Did  the  formal  expression  obtained  for  the  refractive  indices  of 
waves  travelling  along  the  x,  y,  and  z  directions  of  the  orthorhombic  lattice 
guarantee  that  the  optical  behavior  of  the  medium  would  conform  to  the  laws 
of  crystal  optics  and  be  valid  also  for  other  directions  of  light  propagation? 
This  is  answered  in  the  affirmative  sense  in  Part  I,  Section  10. 

(2)  Is  the  birefringence  found,  a  major  or  negligible  effect?  It  was  not 
easy  in  1911  to  carry  out  such  a  comparison  because  in  spite  of  the  existence 
oi  an  elaborate  geometrical  theory  of  crystal  structure  the  question  was 

left  wide  open  what  kind  of  particle  -  molecules  or  groups  .‘f  molecules  - 
was  to  be  placed  in  the  assigned  positions.  The  birefringence  might  well 
have  its  origin  in  an  inherent  anisotropy  of  these  particles  rather  than  in 
the  contribution  caused  by  their  anisotropic  arrangement.  What  I  needed 
for  the  comparison  was  a  crystal  of  known  birefringence  which  was  built 
according  to  the  simple  orthorhombic  lattice  postulated  in  my  calculations. 
The  famous  oJd  professor  of  mineralogy  and  crystallography,  Paul  Groth, 
to  whom  I  took  my  question,  exclaimed  after  some  hard  thinking:  ’There  is 
one  crystal  which  is  most  certainly  built  according  to  the  lattice  you  want, 
for  it  has  three  marked  cleavage  planes  at  right  angles  to  one  another. 

This  crystal  is  anhydrite,  CaSO^.  Try  that.  ’*  So  I  did.  The  result  of  my 
calculations  is  to  be  found  in  Part  I,  Section  11.  It  was  shattering  my  secret 
hopes  of  obtaining  values  showing  at  least  some  similarity  to  the  actual  bire¬ 
fringence  of  anhydrite. 

So,  with  some  resignation,  I  stated  that  the  order  of  magnitude  of  the 
birefringence  due  to  the  lattice  arrangement  is  the  same  as  the  observed 
one.  From  this  it  followed  that  in  any  future  calculation  of  birefringences 
the  effect  of  the  anisotropic  arrangement  would  have  to  be  taken  into  account. 

It  was  only  years  later,  in  1925  and  1926,  that  the  true  structure  of  an¬ 
hydrite  was  established  as  consisting,  essentially,  of  SO^  tetrahedra  ar¬ 
ranged  according  to  a  b-face  centered,  not  a  primitive,  orthorhombic  lat¬ 
tice. 
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In  my  thesis,  I  made  some  attempt  to  discuss  what  would  be  implied 
in  an  extension  of  the  theory  to  more  complicated  assumptions  r«.gu*dfi\g 
the  crystal  structure.  The  two  sections,  one  on  the  influence  of  infrared 
proper  vibrations,  and  the  other  on  the  refraction  by  a  general  periodic 
arrangement  of  dipoles,  did  not  lead  to  any  physical  results  except  showing 
that  some  of  the  conclusions  reached  under  the  simpler  assumptions  would 
be  invalidated:  they  were,  therefore,  omitted  from  the  paper  published 
in  the  Annalen  der  Physik. 

The  calculation  of  double  refraction  became  much  more  realistic  once 
the  determination  of  crystal  structures  had  become  possible  after  Laue’s 
and  Bragg's  work  in  1912.  In  particular,  E.A.  Hyileraas  was  very  success¬ 
ful  in  accounting  for  the  birefringence  of  the  tetragonal  calomel  (Hg0Cl2)4: 
and  of  the  trigonal  Quartz  \SiO£,**  by  an  extension  of  the  theory  to  more 
general  structures  than  the  simple  dipole  lattice.  In  the  case  of  quartz, 
Hyileraas  could  also  account  for  the  optical  rotatory  power,  using  an  ex¬ 
tension  of  my  theory  which  had  meanwhile  been  worked  out  by  M.  Born  in 
order  to  derive  optical  activity  from  the  presence  of  screw  axes  in  the 
crystal  structures.  In  all  these  cases,  it.  was  sufficient  to  assume  the  atoms 
to  be  represented  by  isotropic  dipoles;  later,  cases  became  known  where  an 
inherent  anisotropy  of  the  resonators  had  to  be  assumed  in  addition  to  the 
effect  of  their  anisotropic  arrangement.*** 

The  theory  of  'Idspersion  and  Double  Refraction  of  Electron  Lattices 
(Crystals),"  was  first  printed  under  this  title  in  Gottingen  1912  and  sub¬ 
mitted  as  Ph.  B.  thesis  to  the  Philosophical  Faculty  (IL  Section)  of  the 
University  of  Munich  under  the  date  of  16th  February  1912,  The  present 

*  Zs.  f.  Physik  1926,  36:859-896. 

.  ■  ..  — -  -  —  MA 

**  Zs.  t  Phys.  1927,  14:871. 

-  -  AM 

***  See,  for  instance,  H.  C.  Bolton,  W.  Fawcett,  I.  D.  C.  Garney, 

Proe.  Phys.  Soc.  (1962)  80:199-208  and  (1963)  82:33-46. 

-  'V*W\ 
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translation  is  from  a  subsequent  publication*  where,  under  a  new  general 
title  Foundations  of  Crystal  Optics,  an  introduction  is  given  not  only  to  the 
theory  of  dispersion,  but  also  to  Part  n  which  is  called  Theory  of  Reflec-  - 
tion  and  Refraction.  Two  more  papers  appeared  under  the  same  general 
title,  namely  in  1917**,  and  twenty  years  later.  *** 

The  latter  two  parts  deal  with  the  case  that  the  wavelength  of  the 
optical  field  is  shorter  than  the  distance  between  neighboring  atoms  in  the 
crystal.  The  strictly  valid  transformations  of  the  sum  of  spherical  dipole 
waves  from  Part  I  are  again  used,  but  their  discussion  proceeds  along 
entirely  different  lines  for  the  X-ray  than  for  the  optical  case.  The  same 
condition  of  consistency  of  the  dipole  oscillations  with  the  field  created  by 
them  is  applied  and  its  result  is  the  determination  of  the  mode  of  propaga¬ 
tion  of  the  optical  field  in  the  interior  of  an  unbounded  crystal.  The  X- 
opiical  field,  however,  is  more  complicated  than  that  of  visible  light,  for 
it  consists  not  of  a  single  plane  wave,  but  of  a  pencil  of  such  waves  which 
are  coupled  together.  The  dynamical  condition  for  the  propagation  oi  the 
X-optical  field  is  expressed  geometrically  by  a  surface  of  a  similar  nature 
to  the  surface  of  normals  in  visible  optics  (Part  I,  Section  10) ;  for  X-rays 
this  is  called  the  surface  of  dispersion.  Once  this  surface  has  been  deter¬ 
mined,  all  possible  modes  of  propagation  of  self-consistent  X-fields  in  the 
interior  of  the  unbounded  crystal  are  known.  There  remains,  however,  the 
same  problem  as  in  the  case  of  visible  light,  namely  of  relating  the  internal 
optical  field  to  the  external  one  for  a  crystal  with  boundary,  i.  e. ,  a  crystal 
filling  only  one  half-space.  This  problem  was  treated  for  visible  light  in 


*  Annalen  der  Physik  (Lpz)  (1916)  49:1-38. 

- — -  ,VV\ 

**  Part  HI:  The  Crystal  Optics  of  X-Rays,  Ann.  d.  Phys.  Lpz  (1917) 
54:519-597. 
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♦♦♦Part  IV:  Establishment  of  a  General  Condition  of  Dispersion  in 
Particular  for  X-Ray  Fields,  Zeitsch,  f.  Kristallographie  (A)  (1937) 
97:1-27. 
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Part  II  of  the  series,  and  the  same  principle,  established  there,  could 
later  be  applied  to  answer  the  X-ray  case. 


The  so-called  Dynamical  Theory  of  X-ray  Diffraction  of  Part  IH  was 
restricted  by  the  assumption  that  the  crystal  could  be  represented  by  a 
simple  {,v>rvmitive'>  lattice  of  dipoles.  The  restriction  to  an  orthorhombic 
lattice  was  lifted  because  one  bad  meanwhile  learned  to  handle  ncnorthogonal 
lattices  by  means  of  the  concept  of  reciprocal  vector  sets.  *  hi  the  simple 
case  of  only  two  plane  waves  forming  the  X-opiieai  field,  or. :  could  deal 
with  a  general  crystal  structure,  thus  getting  rid  of  the  fir  a  restriction  {on 
the  crystal)  by  introducing  a  new  one  (on  the  field).  The  i  .‘suiting  intensity 
formula  (hut  not  the  entire  theory)  was  published  in  a  paper  of  mine.  ** 

This,  however,  left  the  general  case  unanswered  in  which  an  X-opticc  1 
field  consisting  of  n  component  plane  waves  travels  through  a  general  peri¬ 
odic  crystal  structure.  This  problem  was  tackled  in  .  art  IV  and  brought 
to  a  general,  if  formal,  solution.  The  ideal  aim  of  ties  theory  was  to 
obtain  a  uniform  treatment  valid  for  short  X-wavei  mgtfcs  as  well  as  for 
the  long  ones  of  visible  light  and  beyond  this  for  the  static  dielectric 
constant,  and  all  that  in  a  quite  general  periodic  medium.  The  formal 
answer  obtained  in  this  paper  deals  only  with  the  internal  condition  of 
self-consistency  which,  however,  is  the  clue  to  the  boundary  problem.  An 
application  of  this  general  theory  to  the  propagation  of  an  X-field  consisting 
of  three  plane  waves  in  a  diamond-type  structure  is  to  be  found  in  the  paper 
by  Ewald  and  Heno.  *** 

The  dynamical  theory  of  X-ray  diffraction  was  recast  in  a  different 
form  in  1931  bv  M.  von  Laue.  The  crystal  was  considered  to  be  a  con¬ 
tinuous  periodic  medium  in  which  a  periodic  dielectric  constant  affects  the 
propagation  of  the  X-rays.  Laue's  theory  thus  departs  from  the  strict 

*  Introduced  by  Willard  Gibbs  in  his  lectures  1881-1884,  but  first 
published  in  Glbbs-Wiison  Vector  Analysis,  1901. 


Physik.  Zeitsch  (1925)  26:29-32,  and  1926,  27:182. 


***  Acta  Crystallographica  1S68,  24:5-15. 


A  9 


concepts  of  the  Lorentz  electronic  theory  of  matter  and  rashes  use  of  the 
Maxwellian  material  constants  and  field  quantities  and  their  relations.  A 
justification  for  this  step  was  given  by  M.  Kohler  in  1935  on  the  basis  of 
wave  mechanics  and  perturbation  theory.  The  chief  difference  between* 

Laue's  and  my  own  version  of  the  theory  is  that  Laue  makes  use  of  Max¬ 
wellian  type  boundary  conditions  at  the  surface  of  a  half-crystal  in  order 
to  connect  the  external  and  the  internal  field,  whereas  in  my  theory  this  is 
achieved  by  the  condition  of  self-consistency  which  has  to  be  satisfied  at  each 
dipole  throughout  the  depth  of  the  crystal.  The  results  are  the  same  in  both 
versions. 

The  dynamical  theory  of  X-ray  diffraction  has  achieved  increasing  importance 
where  crystals  are  used  as  optical  instruments  of  considerable  perfection. 

Of  greater  practical  importance,  however,  is  the  application  of  the  concepts 
created  in  this  theory  to  the  case  of  electron  diffraction.  The  scattering 
effect  of  matter  is  much  stronger  for  electrons  than  for  X-rays,  and 
therefore  the  multiple  scattering  effects  are  more  prominent  with  ^cctruns. 
Since  these  are  included  in  the  ’’dynamical”  theory,  this  theory  has  to  be 
applied  to  electron  diffraction  at  an  earlier  stage  of  refinement  than  is 
necessary  in  X-ray  diffraction.  The  theory  of  electron  diffraction  was 
first  given  by  EL.  A.  Bethe  in  1928,  and  it  uses  many  of  the  concepts  in¬ 
troduced  for  X-r»ys.  It  is  a  scalar  theory  of  diffraction  dealing  with 
Schrodinger's  psi  function  as  the  field  quantity  rather  than  with  the  vectorial 
field  strengths  of  the  optical  theory.  The  theory  is  an  unfailing  guide  in  the 
discussion  of  the  highly  complex  observations  of  electron  diffraction. 

With  the  advent  of  neutron  diffraction,  the  concepts  and  methods  of  the 
dynamical  theory  found  an  important  new  field  of  application.  * 

After  this  short  survey  of  the  wide  field  of  developments  beginnir  g  with 
the  two  parts  of  the  series  which  are  to  be  found  here  in  translation,  let  us 
return  to  Part  n  with  some  comments,  ft  contains  the  connection  bei  ween 
the  Internal  and  the  external  field  in  a  "half-crystal,  ”  i.e. ,  a  lattice  lifting 

*  See  for  instance  C.G.  Shull  (1968)  Pbys.  Rev.  Letters,  21:158  5- 
!589.  - 
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one  halt  space  and  cut  oit  at  the  surface  a  ^  0.  The  half-crystal  is  still 
infinite,  but  it  is  bomrded.  Strictly  speaking,  it  has  no  boundary  is  tie 
sense  of  the  Maxwell ian  theory  for  a  continuous  medium,  and  therefore 
the  theory  does  net  contain  boimdary  conditions  which  would  be  valid  over 
the  entire  plane  z  =  0.  instead,  the  half-crystal  is  an  open  structure  with 
empty  space  not  only  outside  the  lattice,  but  also  inside  between  the  dipoles. 
The  peculiar  boundary  conditions  for  a  continuous  medium,  like  E^^  and 
_  the  same  on  both  sides  of  the  boundary,  axe  replaced  in  the  open 
structure  by  the  condition  of  self-consistency  at  the  site  of  each  dipole,  no 
matter  at  what  depth  under  the  surface  it  lies. 

It  was  already  pointed  out  in  Part  I,  that  in  the  theory  of  dispersion 
there  was  no  room  for  an  incident  wave  -  in  fact  where  should  it  come  from 


if  the  crystal  fills  all  space  ?  In  Part  H,  it  is  shown  that  by  restricting  the 
summation  of  spherical  dipole  wavelets  to  source^  on  one  side  of  the  plane 
z  =  0,  extra  fields  are  created  which  travel  on  both  sides  away  from  this 
plate  with  the  velocity  c,  or  in  the  higher  terms  of  the  series,  are  in¬ 
homogeneous  surface  waves;  i.e. ,  waves  whose  (small)  phase  velocity  is 
directed  along  the  surface  while  their  amplitude  decreases  exponentially 
with  increasing  distance  from  the  surface.  The  term  ’boundary  waves"  is 
used  for  these  waves  in  order  to  encompass  both  the  homogeneous  and  in¬ 
homogeneous  ones,  hi  other  places*  I  have  called  them  epiwaves,  in 
contrast  to  the  field  in  the  unbounded  crystal  which  consists  of  mesowaves. 

In  Part  n,  the  boundary  waves  arise  from  residue  formation  at  special 
series  of  poles  in  the  plane  of  integration.  A  less  formal  understanding  of 
their  occurrence  is  desirable  and  can  be  obtained  in  two  different  ways. 

The  first  is  a  simple  application  of  the  sum  of  a  geometrical  series, 

2 

say  of  quotient  Q.  The  infinite  series,  1  +  Q  +  Q  +  . . .  has  the  sum 
l/(l-Q) ;  the  finite  series,  1  +  Q  +  Q2  +  . . .  has  the  sum  (l-Q^/ 

(1  -Q),  that  is,  an  additional  term  is  caused  by  the  termination  of  the  series. 


*  For  instance  iu  Gomer  and  Smith,  Structure  and  Properties  of  Solid 
Surfaces,  Uolv.  of  Chicago  Press  1953,  Ch.  IL 
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N ow  we  consider  tbs  spherical  waves  issuing  from  a  plane  of  equal  dipole 
phase.  These  waves  will,  soon  after  leaving  that  plane,  merge  to  form  a 
plane  wave  travelling  with  the  velocity  c  away  from  the  plane  of  the  dipoles. 
At  a  field  point  z  there  arrive  such  plane  waves  emitted  from  the  nearest, 
the  next  nearest  anc  successively  more  distant  planes  o?  dipoles.  Since  the 
assumed  velocity  of  the  dipole  phase  is  q,  there  w-fii  be  at  z  a  time  lag  6 
between  contributions  from  neighboring  dipole  planes  equal  to  a(l/q-l/c) 
where  a  is  the  distance  between  the  dipole  planes.  We  thus  obtain  at  z  a 
sum  of  waves  of  equal  amplitudes  but  with  a  phase  factor  Q  =  e 1  between 
them  where  v  is  tbs  frequency  and  )  -  2~  /-T  The  resulting  amplitude  at 
z  is  the  sum  of  the  geometrical  series 

- 

□  n 

If  the  series  is  continued  to  infinite  values  cf  n  (assuming  a  slight  absorption 

for  the  sake  of  convergence),  the  resulting  amplitude  is  1/(1  -Q).  If  we 

st 

have  to  terminate  the  summation  at  the  (N  -  1}  layer  above  z,  an  additional 
term  Q*  /(l -Q)  arises.  This  term  corresponds  to  the  boundary  wave.  A 
full  discussion  of  the  theory  along  these  simplified  lines  is  the  subject  of  a 
paper  of  mine. 

The  second,  even  more  physical  way  of  understanding  the  field  in  the 
interior  of  the  half-crystal  is  this.  Imagine  first  all  space  filled  with  the 
dipole  lattice  oscillating  as  assumed  in  Part  I.  The  result  is  the  field 
derivable  from  the  potential  Part  I,  Section  6,  Eq,  (30)  which  we  new  call 
for  short  the  mesofield.  Next  cancel  the  contributions  of  the  dipoles  above 
the  plane  z  =  0  by  adding  in  the  upper  half  of  space  fictitious  dipoles  of  the 
opposite  phase.  Then  the  total  field  in  the  lower  half-space  is  the  sum  of 
the  mesofield  and  the  field  produced  in  the  lower  half-space  by  the  fictitious 


*  Forts chritte  der  Cherde,  Physik  und  physikal,  Chemie  Serie  B  (1925) 
18,  Heft  8  (Berlin,  Gebr.  Borntrager  Verlag).  See  also  Physica  (1924) 
4:234-251. 
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dipoles,  i.  e. ,  outside  the  fictitious  medium.  This  is,  of  course,  a  field 
consisting  of  waves  of  velocity  c;  in  fact,  it  is  the  field  of  the  internal 
boundary  waves  of  the  half  crystal.  TH  fictitious  dipoles,  apart  from 
destroying  the  mesofield  in  the  upper  half  of  space,  also  create  there  an 
internal  boundary  wave  in  their  (fictitious)  medium.  For  the  real  half- 
crystal  this  is  its  external  boundary  wave  which  then  becomes  the  reflected 
wave  when  the  relationship  to  the  wave  incident  on  the  half-crystal  is 
established. 

This  physical  explanation  is  not  too  different  from  the  foregoing  ex¬ 
planation  by  means  of  the  termination  term  of  the  finite  geometrical  series. 
For  the  formula  for  the  finite  sum  can  he  interpreted  as  the  instruction: 
sum  first  the  infinite  series  and  then  subtract  all  terms  beyond  the  point  of 
termination  (p  ~  1);  i.  e. ,  the  same  infinite  series  but  with  all  terms  multi¬ 
plied  by  a  factor  Qp. 


My  entire  series  of  papers  deals  with  the  optical  problems  strictly  in 

the  sense  of  A.  H.  Lorentz’  Theory  of  Electrons,  but  differs  from  his 

treatment  by  the  fact  that  r.o  averaging  over  random  positions  of  dipoles 

becomes  necessary.  This  is  the  virtue  of  SornmerMd’s  clear  formulation 

of  the  problem.  Shortly  after  ray  work,  and  independent  of  it.  Professor 

C.  W.  Oseen  in  Oslo  found  the  same  screening  action  of  a  boundary  that  is 

so  essential  in  my  treatment.  His  work  is  based  on  Maxwell’s  equations 

for  a  homogeneous  medium.  He  assumes  a  light  wf’-'e  to  fail  on  a  finite 

body  bounded  by  a  surface  S.  on  which  the  usual  boundary  conditions  hold. 

By  transforming  the  integral  expressing  the  polarization  created  in  tbs  body, 

he  shows  that  this  consists  of  the  two  parts,  one  of  which  conforms  to  the 

wave  equation  4u  +  k6  u  =  0  with  x  =  K  =.  uk  and  the  other  with  x  =  k  , 

O  0  o 

where  kQ  =  wave  constant  in  free  space,  n  ~  refractive  index  and  Kq  =  wave 
constant  in  the  medium.  The  same  holds  for  the  total  field,  which  is  the 
field  created  by  the  volume  polarization  in  the  body  plus  the  field  incident 
from  outside.  The  condition  of  self-consistency  in  the  body  can  be  ful¬ 
filled  only  if  those  parts  of  the  field  vanish  that  progress  with  velocity  c.- 
That  is,  the  external  field  into  which  the  body  is  placed,  must  cancel  that 
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part  of  the  field  which  originates  in  the  volume  polarization  and  progresses 
with  velocity  c.  i  xie  latter  part  can  be  expressed  as  a  surface  integral  of 
the  polarisation  over  the  boundary  surface  S.  Thus  it  is  by  the  action  of 
the  surface  that,  the  body  is  shielded  from  the  impinging  optical  field. 

This  is  in  general  terms  the  same  result  that  I  obtained  for  the  special 
case  of  a  plane  boundary.  The  Oseen  theory  is  well  rendered  elsewhere.  * 

Looking  back  on  the  mathematical  technique  used  in  Parts  I  and  n,  it 
may  well  be  considered  cumbersome.  When  the  papers  were  written,  the 
technique  of  Fourier  transformation  was  still  in  its  infancy.  Reciprocal 
sets  of  axes  had  yet  tc  be  discovered  by  physicists.  Their  importance  for 
the  Fourier  development  of  functions  in  a  non-orthogonal  cell  was  shown  in 
my  paper  on  the  reciprocal  lattice,  but  even  then  the  idea  of  considering 
the  Fourier  coefficients  as  the  transform  of  a  periodic  function  had  either 
not  been  pronounced  or  not  gained  the  attention  of  the  physicists.  Now¬ 
adays,  the  transformation  of  the  sum  of  spherical  dipole  waves  into  a  sum 
of  plane  waves  is  achieved  easily  by  applying  a  few  general  properties  of 
Fourier  transformation.  This  method  can  be  found  in  my  paper,  *** 
recording  a  lecture  given  in  1932. 

My  papers  contain,  for  the  first  time  I  believe,  a  clear  separation  of 
the  internal  problem  of  finding  the  kind  of  optical  field  that  can  exist  inside 
a  dispersive  medium,  from  the  external  problem  of  finding  out  what  happens 
if  an  external  wave  falls  on  the  surface  of  such  a  medium.  It  is  only  lately 
that  I  have  become  aware  that  the  splitting  of  the  entire  problem  into  these 
two  parts  corresponds  quite  closely  to  the  dbnsion  into  two  steps  of  the 
mechanical  problem  of  oscillations  generated  by  an  impact  on  a  mechanical 
system.  There,  a  preliminary  study  is  made  of  the  free  oscillations  of  the 
system,  and  the  knowledge  sc  gained  is  subsequently  applied  to  the  behavior 

*  Rom  and  Wolf,  Principles  of  Optics,  2nd  ed.  (1964)  pp.  100-104, 
Pergamon  Press.  The  original  paper  by  Oseen  is  Ann.  d.  Phys.  (Lpz.), 
(1915)  48mp.  1-56. 

**Zeitsch,  f.  Kristallog  (1921)  56:129-156, 

. .  AAA 

***Anoales  de  l’institut  Henri  Poincare,  8:79-110. 
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of  the  system  after  an  impact  has  set  it  in  motion.  The  free  oscillations 

are  those  requiring  no  outside  forces  and  they  can  exist  ccly  at  certain 

frequencies.  The  theory  of  dispersion  (Part  I;  omits  all  external  forces; 

to  the  determination  of  the  proper  frequency  (or  frequencies)  of  the 

mechanical  system  there  corresponds  the  determination  of  the  wave  vector 

K  (or,  in  the  case  of  X-ravs  of  the  several  wave  vectors  K,^)  for  which 

a 

alone  the  system  of  field  and  dipoles  is  self-consistent.  This  knowledge  is 
used  to  construct  the  state  of  the  system  filling  a  half-space  when  an 
impact  is  given  to  it  by  an  externa]  wave  impinging  on  its  surface.  The 
development  of  the  mechanical  system  in  time  is  paralleled  by  the  spatial 
development  in  the  optical  system.  This  analogy  will  be  found  worked  out, 
for  the  case  of  X-rays,  in  a  recent  paper  of  mine.  * 

It  is  deplorable  that  most  physicists  interested  in  optics,  as  a  sub¬ 
division  of  Physics,  ignore  or  nearly  ignore  the  fascinating  branch  of 
X-ray  optics.  In  my  papers,  visible  and  X-ray  optics  are  joined  together 
in  a  comprehensive  theory  which,  because  of  the  mathematical  strictness 
of  its  foundation,  holds  for  the  whole  range  of  wavelengths,  from  the  static 
case  (A  Infinite),  through  the  visible  (A  ~  5000  to  X-rays  (A~1A)  . 

The  discussion,  it  is  true,  varies  in  these  different  regions,  but  the 
principles  on  which  the  theory  :s  based  remain  the  same.  This  wide  range 
has  attracted  me  to  the  theory  through  about  SO  years-  in  the  course  of 
which  the  X-ray  optics  gradually  developed  from  a  purely  logical  necessity 
to  a  reliable  guide  through  a  wealth  of  experimental  observations  which 
became  feasible  with  the  advent  of  the  art  of  growing  nearly  faultless 
crystals. 


*  Acta  Cryst.  (1969)  A25:103-!08. 
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Appendix  B 

Biography  of  P.P.  Ewoid 


Paul  P.  Ewaid  was  born  23  January  138k  in  Berlin,  Germany.  His 
faUier,  a  promising  young  historian,  died  of  appendicitis  at  the  age  of  37 
shortly  before  he  was  bore.  His  mother,  being  free  to  travel,  took  him  on 
extended  trips  tc  England  and  to  Paris.  Thus  by  the  age  of  five  he  had  a 
rudimentary  knowledge  of  both  English  and  French. 

He  graduated  from  the  Victoria  Gymnasium  of  Potsdam  obtaining  his 
Arbitnr  (school-leaving  examination)  in  the  fail  of  1805.  Tee  vsm  winter 
terms  1905/06  were  spent  at  Gonvxlle  and  Caius  College,  part  of  the 
University  oi  Cambridge,  England;  1906/07  at  the  University  of.  Gottingen. 
From  1907  to  1912,  he  studied  under  Professor  Arnold  Somme  rie  Id  at  the 
University  of  Munich  and  received  his  Ph.  D.  in  Physics  in  February  1812. 
ms  choice  of  a  thesis  subject  influenced  Ms  entire  scientific  career, 
^Dispersion  and  Double  Reiraetlon  of  Electron  Lattices  ( Crystal  o)"  dealt 
with  the  propagation  of  light  in  crystals  and  was  explored  by  him  further  in 
various  ways  in  later  papers. 

The  years  of  World  War  I  were  spent  as  an  X-ray  technician  of  the 
German  army  on  the  Russian  front. 

Returning  to  university  life  at  the  University  of  Munich  in  1918,  he 
became  Privatdozent  (lecturer)  there,  but  soon  moved  to  the  Technische 
Hoctschule,  Stuttgart,  in  1921  as  Professor  of  Theoretical  Physics,  a 
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position  he  held  until  1937.  During  this  time  he  wrote  a  book:  Kristalle 
und  Rontgenstrahlen  (Crystals  and  X-rays),  (Springer  1923),  and  together 
with  C.  Hermann,  Strukturbericht  1913-1928  (Structure  Report  1913-1928), 
(Akad,  Verlagsges.  1931). 

His  first  visit  to  the  United  States  occurred  in  1936  when  he  was 
invited  to  give  a  course  of  lectures  on  X-ray  diffraction  at  the  Summer 
School  in  Physics  in  Ann  Arbor.  Leaving  Germany  in  1937,  he  returned 
to  Cambridge  University,  England,  where  from  1937  to  1939  he  was  a 
Research  Fellow  in  the  Crystallographic  Laboratory.  Then  he  moved  to 
The  Queen's  University,  Belfast,  first  as  Lecturer,  then  as  Professor  of 
Mathematical  Physics. 

After  ten  years  in  Belfast,  he  moved  to  the  United  States  and  in  1949 
became  Professor  of  Physics  and  Head  of  the  Physics  Department  at 
Polytechnic  Institute  of  Brooklyn,  Brooklyn,  New  York.  Jr.  1957  he  re¬ 
linquished  the  position  as  head  of  the  physics  department  but  remained  as 
Profeseor  of  Physics  until  he  retired  in  1959. 

Over  his  long  and  productive  career  he  has  written  numerous  papers 
of  which  he  considers  some  to  be  fundamental.  He  is  known  for  his  Dart  in 
the  Ewaid-Oseen  extinction  theorem,  and  Ewald's  reciprocal  lattice  is  a 
subject  of  study  for  students  of  optics  and  crystals.  It  is  an  important 
tool  for  all  work  on  X-ray  diffraction  in  crystals.  Ewald's  principal 
achievement  is  the  "Dynamical  Theory  of  X-ray  Diffraction”  which  he  dealt 
with  in  a  number  cf  papers  using  the  same  principles  he  lad  employed  in 
Ms  doctoral  thesis  and  his  two-part  paper  published  in  1916  of  which  his 
thesis  was  Part  L  These  papers  elucidate  the  principles  which  govern  the 
transmission  of  electromagnetic  radiation  through  a  crystalline  medium. 
The  wavelength  range  which  is  considered  includes  both  visible  light  and 
X-rays.  Because  of  his  work  he  has  been  the  recipient  of  a  number  cf 
honors.  He  was  awarded  an  honorary  doctorate  from  the  Technische 
Kcchschule,  Stuttgart  in  1954;  University  of  Paris  in  1958;  Adelpni  Univer¬ 
sity,  New  York  in  1966;  and  the  University  of  Munich  ia  1967.  He  was 
made  a  Fellow  of  the  Institute  of  Physics  and  Physical  Society  (England) 


in  1945;  the  American  Physical  Society,  13£2;  the  American  Academy  of 
Arts  and  Sciences,  1955;  The  Royal  Society,  1958;  Akademic  Deutsche? 
Naturforscher  (Leopoldina)  1966. 

He  is  corresponding  member  of  the  Akademic  oer  Wissensebaffen  zu 
Gottingen  (1937)  and  of  the  Bayerische  Akademie  der  Wissenechaffen  in 
ftfiincfaen,  (1962).  He  was  made  an  honorary  iremoer  of  Societe  fran9aise 
de  mineralogie  et  cristallographie  in  1955;  Deutsche  Mineraloglsche 
Gesellschaft  1958;  Cambridge  Philosophical  Society,  1969. 

Over  the  years,  in  addition  to  his  numerous  papers,  he  found  time  to 
be  a  co-editor  of  Zeitschr.  f.  Kristallographie  from  1922  to  1937;  and  to 
be  editor  of  Acta  Crystallographies,  from  its  earliest  beginnings  in  1948  to 
1959.  In  addition  to  the  books  previously  mentioned,  be  contributed  sec¬ 
tions  to:  Landolt- Borns tein  Tabellen  (two  editions);  Muller-Pouillet  Lehr- 
buch  d.  Physik  (11th  Edition  1926)  (An  English  translation  is  in  'The  Physics 
of  Solids  and  Fluids"  (1st  Ed.  1930,  2nd  Ed.  1936)  (Blackie));  Hasdbuoh  d« 
Physik  (Springer  1927  and  1933) ;  High  Speed  Aerodynamics  Vol.  1  (Editor, 

F.  D.  Rossini,  Princeton  University  Press  1955).  He  was  the  principal 
author  as  well  as  editor  of  the  book  "Fifty  Years  of  X-ray  Diffraction” 
published  by  the  International  Union  of  Crystallography  in  July  1962  to 
celebrate  the  Laue  and  Bragg  discoveries. 

He  was  Member  of  the  executive  Committee  of  the  International  Union 
of  Crystallography  from  1948  to  1988,  and  President  of  the  same  organiza¬ 
tion  from  1960  to  1963. 

The  January  1968  issue  of  Acta  Crystallographies  was  dedicated  to  him 
on  the  occasion  of  his  80th  birthday.  The  issue  contains  one  of  his  scientific 
articles  and  a  three  page  article  with  the  title  ’Personal  Reminiscences", 
which  contains  many  interesting  biographic  remembrances. 

His  retirement  from  university  life  in  1959  was  not  a  retirement  from 
other  activities.  True,  he  moved  to  quieter  surroundings  in  Xew  Milford, 
Connecticut,  but  as  indicated  by  dates  already  given  he  has  remained  active. 
For  a  number  of  years  he  maintained  a  rather  busy  schedule  of  lectures  in 
addition  to  his  activities  for  the  International  Union  of  Crystallography. 


Presently,  (May  1970),  he  continues  an  active  interest  in  science  and  car¬ 
ries  on  a  scientific  correspondence  of  some  magnitude  and  occasionally 
travels  to  a  meeting  even  if  it  is  far  away  from  Connecticut. 
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